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This  report  develops  a  technique  for  bounding  the  maximum  voltages 
and  currents  at  terminations  of  a  multiconductor  transmission  line  (MTL) 
located  behind  an  aperture-perforated  conducting  screen  excited  by  an 
electromagnetic  pulse.  The  electromagnetic  field  Is  coupled  through  a 
small  aperture  to  provide  the  excitation  of  a  multiconductor  transmission 
line  behind  the  aperture.  A  model  is  presented  in  terms  of  external  and 
internal  sources  which  in  tun  creates  traveling  waves  on  the  mulci- 
conducuor  transmission  line.  The  latter  transfer  energy  to  the  termina¬ 
tions.  The  energy  at  a  termination  is  translated  to  voltages  and  cur¬ 
rents  from  which  the  upper  bounds  are  determined.  These  upper  bounds 
are  obtained  using  vector  norms  and  associated  matrix  norms.  The 
formulation  is  presented  in  the  frequency  domain  and  transformed  to 
the  time  domain  to  obtain  useful  upper  bounds  for  transient  analysis 
of  multiconductor  transmission  line  geometries  with  aperture  excitation. 
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CHAPTER  I 


INTRODUCTION 

In  designing  some  system s,  Che  designer  should  be  able  to  charact¬ 
erize  the  penetration  of  electromagnetic  pulses  (EMP)  or  lightning 
signals  through  apertures  of  general  shapes  as  well  as  quantify  the 
effects  of  the  coupled  energy  on  transmission  lines  located  in  the 
vicinity  of  the  aperture. 

Apertures  that  are  of  concern  to  the  designer  are  usually  electro- 
magnetically  small  over  the  spectrum  of  the  EMP,  or  lightning,  and  their 
existence  may  be  for  some  purpose,  e.g.,  windows,  open  access  holes,  or 
they  may  be  unintentional  as  in  the  case  of  cracks  around  doors  or 
plates  covering  access  ports  or  poor  electrical  seams  [1].  Small  in 
the  sense  of  electromagnetic  penetration  implies  that  the  maximum 
dimension  of  Che  aperture  is  small  compared  with  the  wavelength  of  the 
time-harmonic  electromagnetic  field.  The  analysis  of  the  coupling 
(penetration)  problem  has  Been  investigated  by  a  large  number  of  people 
since  1897.  The  first  scientist  to  propose  a  solution  was  Lord 
Rayleigh  [2,3],  whose  solution  was  expressed  as  an  ascending  power 

2n 

series  of  the  wavenumber  k(.“— )  where  X  is  the  wavelength.  Bethe  [4] 
presented  the  results  for  the  leading  terms  in  the  Rayleigh  method  by 
using  a  scalar  potential  function  approach.  Bouwkamp  [5]  used  a  set  of 
coupled,  integro-dif ferential  equations  to  solve  the  problem.  Recently, 
some  have  used  an  integral  equation  approach  to  tackle  the  problem. 
Notable  among  them  are  3utler  and  Umashankar  [6], [7]. 
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The  coupling  of  che  energy  from  an  incident  electromagnetic  wave 
to  a  transmission  line  located  behind  an  aperture-perforated  conducting 
screen  has  been  investigated  by  many  engineers  and  physicists  in  the 
past  decade.  Kajfez  [8]  has  computed  the  coupled  energy  by  the  use  of 
equivalent  electric  and  magnetic  dipole  moments .  He  has  derived  the 
equivalent  sources  of  a  transmission  line  model  by  using  both  mode¬ 
matching  and  reciprocity  techniques. 

Butler  and  Umashankar  1 9 J  have  approached  the  problem  numerically 
by  the  method  of  moments,  and  have  formulated  integro-differential 
equations  for  a  finite-length  wire  with  arbitrary  orientation  behind  an 
arbitrarily  shaped  aperture.  They  have  also  taken  into  consideration 
the  scattered  energy  into  the  aperture  by  the  wire  after  wire 
excitation.  Lee  and  Tang  110]  have  solved  the  problem  using  transform 
techniques,  and  have  obtained  the  same  equivalent  sources  as  Kajfez  [8]. 
In  addition  they  have  determined  the  effects  of  a  wire  being  very  close 
to  the  aperture.  Davis  [11]  has  developed  a  model  using  a  method 
spatially  equivalent  to  Lee  and  Yang  [10]  method.  He  has  also  come  up 
with  a  capacitive  term  in  his  model  which  is  not  in  Lee  and  Yang.  The 
importance  of  che  capacitance  occurs  when  the  wire  is  close  to  the 
aperture  causing  a  capacitive  discharge  of  the  aperture  region  by 
currents  on  the  wire. 

Davis  [11]  has  also  found  a  method  for  bounding  the  maximum  voltage 
and  current  levels  at  terminations  of  a  wire  behind  an  aperture.  In  his 
analysis,  he  has  approached  the  problem  both  in  the  frequency  and  time 


domain. 
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This  report  extends  the  bounding  problem  of  a  single  wire  to  the 
problem  of  obtaining  an  upper  bound  for  the  computation  of  the  voltages 
and  currents  at  terminations  of  multiconductor  transmission  lines  (MTL) 
located  behind  an  aperture-perforated  conducting  screen.  The  computa¬ 
tions  are  carried  out  in  both  frequency  and  time  domains.  There  are 
several  stages  that  lead  to  the  final  results. 

In  Chapter  II,  the  electromagnetic  field  coupling  through  small 
apertures  is  discussed.  This  is  then  extended  in  Chapter  III  to  the 
excitation  of  a  MTL  where  a  model  is  presented  in  terms  of  external  and 
internal  sources  which  in  turn  create  traveling  waves  on  the  MTL.  These 
waves  transfer  the  coupled  energy  to  the  terminations.  The  available 
energy  found  at  the  terminations  is  translated  to  voltages  and  currents 
by  impedance  transformations.  Finally,  the  bounds  on  the  voltages  and 
currents  are  obtained  in  Chapters  IV  and  V  using  vector  norms  and 
associated  matrix  norms.  The  discussion  is  closed  with  the  presentation 
of  comprehensive  examples  which  elaborate:  on  the  use  of  the  techniques 
to  find  upperbounds  on  the  signal  levels  in  both  the  frequency  and  time 


domain. 


CHAPTER  II 


ELECTROMAGNETIC  FIELD  COUPLING  THROUGH  SMALL 
APERTURE  IN  A  CONDUCTING  SCREEN 

In  this  chapter  the  equivalent  electric  and  magnetic  dipoles 
representation  of  an  aperture  is  illustrated  and  the  basic  formulas 
are  given  for  this  problem.  The  presentation  is  tutorial  and  the 
scope  is  limited  to  the  discussion  of  techniques  which  may  be  en¬ 
countered  in  application.  The  development  follows  the  procedures 
of  3utler  [13 J. 

Ideally,  the  problem  will  be  discussed  for  a  small  aperture  in 
an  infinite  planar  screen.  Of  course,  in  practical  problems,  one 
never  can  encounter  an  infinite  or  totally  flat  screen.  However,  if 
the  following  conditions  are  satisfied,  then  one  can  replace  the 
real  world  problem  with  an  ideal  one: 

1.  The  minimum  distance  across  the  surface  of  the  screen 
should  be  large  relative  to  the  wavelength  in  the  medium. 

2.  The  minimum  radius  of  the  curvature  of  the  screen  should 
be  large  compared  with  the  wavelenght. 

3.  The  point  of  observation  should  be  close  to  the  surface 
of  the  screen  relative  to  the  minimum  distance  across  the 
screen. 

4.  The  point  of  observation  should  not  be  close  to  the  outer 
edge  of  the  screen. 

5.  The  edges  of  the  aperture  should  not  be  close  to  the  edges 
of  the  screen. 
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In  general,  the  problem  of  concern  is  shown  in  Figure  (1),  where 
one  sees  an  infinite  conducting  screen  with  an  aperture  (A)  cut  in 
it.  The  screen  is  extremely  chin  and  separates  the  space  into  two 
parts,  each  characterized  by  u  (medium  permeability)  and  e  (medium 
permittivity)  where  e  can  be  complex  to  account  for  a  lossy  medium. 
On  each  side  of  the  screen,  there  are  electric  and  magnetic  sources 
J+  and  M+  which  are  known  impressed  currents,  and  give  rise  to  the 
excitation  of  the  aperture. 

One  can  show  [14]  that  the  diffracted  field  by  a  3mall  aperture 
can  be  represented  approximately  by  radiation  from  equivalent  elec¬ 
tric  and  magnetic  dipoles.  Figures  (2)  and  (3)  depict  the  idea  of 
electric  and  magnetic  dipole  representation  for  an  aperture.  Mote 
that  in  both  cases  the  equivalent  dipoles  radiate  in  the  presence 
of  screen  when  aperture  is  completely  closed.  Also,  observe  that 
the  electric  dipole  is  in  the  same  direction  of  incident  electric 
field,  whereas  the  magnetic  dipole  is  in  opposite  direction  of  the 
incoming  magnetic  field. 

The  moments  of  these  equivalent  electric  and  magnetic  dipoles 
are  related  to  the  known  exciting  fields  by  special  constants  called 
polarizabilities.  (See  Appendix  A  for  details.)  For  the  case  of 
Figure  (1),  one  can  replace  the  aperture  (A)  by  equivalent  dipoles, 
and  the  moments  of  these  dipoles  can  be  related  to  the  polarizabi¬ 
lities  [11]  as 

p"*p*z*  +  ea  ESc-(r  )  -  E^c"r(r  )jz 
e  ©  ■"  ©  ^_  z  o  -  o  J 


(1) 


PEC 


(1) .  Location  of  sources  with  respect  to 

the  aperture  and  the  perfect  electric 
conductor  (PEC) . 


(d) 


(f) 


(•) 


Figure  (2) .  Equivalent  electric  and  magnetic  dipole  representation 
of  an  aperture  [13] :  (a)  electric  field  incident  on 

screen  without  aperture,  (b)  electric  field  incident 
on  screen  with  aperture,  (c)  equivalent  electric 
dipole  on  screen  with  closed  aperture,  (d)  magnetic 
field  parallel  to  screen,  (e)  magnetic  field  and  screen 
with  aperture,  (f)  equivalent  magnetic  dipole  on  screen. 
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Aperture 
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igure  (3).  Equivalent  problem  description:  (a)  incident 
field  on  aperture,  (b)  equivalent  dipole 
moments  replacing  the  aperture. 


where  pg“  is  Che  moment  of  Che  equivalent  electric  dipole,  E^  ~  is 

called  the  short-circuit  electric  field  (field  which  would  exist  if 

there  was  no  aperture) ,  rQ  is  the  point  of  evaluation  located  in 
_cc+  - 

Che  aperture,  and  H  (r  )  is  the  short-circuit  magnetic  field. 

o 

The  plus  and  minus  signs  represent  the  regions  of  the  positive  and 
negative  2  respectively. 

In  addition,  the  constants  a  and  a  are  the  electric  and 

e  m 

magnetic  polarizabilities  respectively;  the  latter  is  a  dyad  defined 

as  a  *  (a  )xx  +  (a  )yy.  Thus  Equation  (2)  can  be  written 
m  m,xx  m,  yy  J 

in  component  from  as 


+ 

1“ 

-  + 

a  ]hSc' 

m,x 

m,xx[_x 

± 

f7.Sc 

r 

■  + 

a  IH 

m,y 

m,yyi_y 

•  + 

+ 

«  + 

1  ~  » 

1  p" 

x  +  p~ 

m 

m. 

x  rm. 

(3a) 

(3b) 

(3c) 


The  equivalent  dipoles  are  located  at  vq  on  the  screen  with  the 

pair  p*  p  + placed  on  the  positive  side  of  the  closed  surface,  and 
e  m 

pg,  p^  placed  on  the  negative  side  of  the  closed  surface. 

As  seen  from  Equation  (1) ,  p  “  is  perpendicular  to  the  screen 

-  ±  -Set  - 

surface  while  the  direction  of  p  is  determined  by  H  (r  )  and 

m  o 

the  shape  of  the  aperture  through  a  .  Now,  at  the  present  stage, 

m 


one  is  able  to  calculate  the  field  which  passes  through  a  hole  in 
a  conducting  surface  as  well  as  the  field  which  scatters  back. 
Consider  Figure  (3)  as  depicted  in  the  two  parts.  Note  that  the 
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direction  of  two  electric  dipoles  are  opposite  to  each  other  in  two 
regions  as  are  the  direction  of  two  magnetic  dipoles. 

Considering  Figure  (3b),  one  can  write: 


E  -(r)  =■  ESc±(r)  +  E6*^)  +  !”*(?)  (4a) 

H  ±(r)  -  HSc~(r)  +  £**(9)  +  Hm^(r)  (4b) 


where  E^  (r)  and  S*(r)  are  total  electric  and  magnetic  fields  present 
at  point  (r)  in  region  II  or  I  corresponding  to  the  (+)  and  (-) 
respectively.  The  (e) and  (m)  superscripts  denote  the  fields  of  the 
electric  and  magnetic  dipoles  in  the  presence  of  the  plane  con¬ 
ductor,  and  the  Sc  denotes  the  short-circuit  fields  defined 
previously. 

The  electric  and  magnetic  fields  at  a  point  (r)  due  to  electric 
and  magnetic  dipoles  located  at  (r^)  in  a  homogeneous  space  of 
infinite  extent  are  given  by  [13,15,16] 


E*  »  -  i  ?x(pexVG) 

(5a) 

H6  ■  -  jo:  PgX^G 

(5b) 

Em  *  ju>  u  p  x7G 
m 

(5c) 

-  -  7x(p  *7G) 
tn 

(5d) 

where  G  is  the  free  space  Green's  Function  defined  as: 
-jk|r  -  rl 


G(r,ro) 


4  ir  I  r  -  r 


e 


(6) 
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Substituting  Equations  (5)  into  (4)  and  using  image  theory,  one 
obtains: 

1  ±G)  -  ESc~  (r)  -  7  7x[p^xVG(r,r)  ] 

+  2ju  p  p  -x7G(r,rJ  (7) 

m  o 

H  ±(r)  -  HSc-(r)  -  2ju>pg-xVG(r  ,rQ) 

-  2Vx [p  “x7G (r , r  ) ]  (8) 

m  o 

If  aperture  reference  is  centered  at  the  coordinate  origin,  then 
rQ  becomes  a  zero  vector  in  the  above  equations.  The  above  equa¬ 
tions  are  valid  in  the  distant  region  of  the  aperture  (at  least  one 
aperture  dimension  away) . 

The  merit  of  the  above  approximations  for  computing  fields  de¬ 
pends  on  the  electrical  size  of  the  aperture,  the  distance  of  observa¬ 
tion  point  from  the  aperture,  and  the  choice  of  the  coordinate 
origin  with  respect  to  which  the  dipole  moments  are  calculated. 

The  details  of  dipole  approximation  are  given  in  Appendix  (B) . 


CHAPTER  III 


ELECTROMAGNETIC  EXCITATION  OF  MULTI CONDUCTOR 
TRANSMISSION  LINES  THROUGH  AN  APERTURE-PERFORATED 
CONDUCTING  SCREEN 


In  chis  chapter,  the  electromagnetic  excitation  of  multiconductor 
transmission  lines  (MTL)  located  behind  an  aperture-perforated  conduct¬ 
ing  plane  is  developed  using  the  procedures  and  method  of  Kajfez  [12]. 
In  order  to  do  this  development  systematically,  an  aperture  representa¬ 
tion  by  dipoles  is  discussed  which  notationally  is  different  from  the 
representation  introduced  in  chapter  II.  The  equivalent  source  models 
are  then  derived  using  methods  introduced  by  Kajfez  [12]. 

Consider  Figure  (4)  with  a  multiconductor  transmission  line 
parallel  to  the  aperture  on  the  plane.  There  exists  an  electric  field 
in  the  aperture  region  which  may  be  designated  by  EA  (z,x).  The 
aperture  may  then  be  replaced  with  the  equivalent  magnetic  surface 
current  density  given  by 


Js  (x,z)  -  (x,z)  x  y 

where  y  is  the  normal  unit  vector  into  the  region  of  interest  (y  >  o), 
For  coupling  to  a  MTL  by  a  small  aperture  this  magnetic  current 
distribution  may  be  replaced  by  the  two  current  dipoles  just  above 
the  closed  aperture  as  shown  in  Figure  (5).  These  two  currant 
dipoles  have  amplitudes  defined  bv 


f r  ' 

' '  Jcm  (x,z)  dxdz 


(9) 


and 


juj£  j  i  (xx  +  zz)  X  Jg  (x,  z)  dxdz. 


(10) 
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The  derivation  of  Equations  (9)  and  (10)  can  be  found  in  Appendix  (B) . 

The  dipole  moments  pg  and  as  introduced  in  Chapter  II  are  also 
proportional  to  the  above  weighted  integrals  of  the  equivalent  surface 
magnetic  current  [11].  From  this  proportionality,  we  may  relate  the 
current  dipole  amplitudes  and  dipole  moments  for  e^uC  by 

*  j«  P0  (11) 

and 

K  a  •  (i2) 


The  concept  of  polarizability  (Chapter  II)  can  be  invoked  here  to 
express  cg  and  c^  in  terms  of  the  exterior  fields.  Considering 
Figures  (6a)  and  6(b),  for  the  case  of  coupling  to  quasi-TEM  (Transverse 
Electromagnetic)  waves,  only  the  y  component  of  E  and  the  x  component 
of  H  are  of  interest.  The  resultant  expressions  are  given  by 


m 


c  x 
mx 


-juuam,xxCHxSC~-  HxSC+  >  *  (13) 


c  y 


Sc-  -  Sc+  .  ~ 
juea  (E  -  E  ) v 
Jey  7  y 


(14) 


Having  introduced  the  preliminary  notation  t  it  is  possible  now  to 
obtain  equivalent  sources  for  a  MTL  behind  an  aperture-perforated 
screen.  Figure  (3)  may  be  redrawn  for  a  different  perspective  as 
shown  in  Figure  (7).  Using  the  concept  of  traveling  waves  introduced 
in  Appendices  (D)  and  (E) ,  the  distribution  of  the  electric  and  magnetic 
fields  for  the  ith  mode  on  the  MTL  traveling  in  (+z)  direction  are 
E1(x,y,z)  ■  ai  e~^SiZ  e^(_x,y)  (15) 

aT(x,y,z)  *  ai  e  ^3lZ  hi(x,y)  .  (16) 


Figure  (6a).  Plane  (p)  with  aperture  (A). 


Figure  (6b).  Aperture  replaced  by  a  tnetalic  lid 
and  the  associated  current  dipoles. 


L-J 

Aperture 


Aperture  replaced  by  a  metalic 
lid  and  associated  current 
dipoles  appropriately  located. 


where  e^  (x,y)  and  h^Gc,y)  are  Che  power  normalized  modal  fields 
over  che  cross  seccion  of  che  MTL.  As  depicted  in  Figure  (7),  ic  will 
be  assumed  chat  chere  are  cwo  inf inices imally  spaced  planes,  denoted 
by  (R)  and  (L) ,  at  each  side  of  che  aperture  junction.  The  equivalent 
source  model  of  MTL  coupling  will  be  established  between  these  planes. 

According  to  Appendix  (F) ,  the  pair  c^  and  c  excite  the  ith 
mode  traveling  in  (+z)  direction  as 

a  .  *  \  [c  h  (x  ,0)  -  c  e  .(x  ,0)]  (17) 

si  2  rnx  xi  0  ey  yi  0 


where  agi  is  the  source  of  a  traveling  mode  in  (+z)  direction.  The 
quantity  hxi(xgj°)  is  the  (x)  component  of  the  ith  modal  magnetic-field 
distribution  evaluated  at  the  point  (x  *  x^)  and  (v  *  0) .  Similarly 
e^fXgjO)  is  the  y-component  of  the  normalized  electric  modal  field. 

A  traveling  wave  source  propagating  in  the  (-z)  direction,  due  to  the 


pair  is 


b  .  *  -  ^[c  h  (xn,0)  +  c  e  . (xo,0)] 
si  2  mx  xi  0  n» 


ey  yi  O’ 


Figure  (8)  shows  a  .  and  b  .  sources  which  are  convenient  for  analvsis 
si  si 

using  scattering  coefficients.  Using  the  Dirac  notation  of  Appendix 


(C) ,  one  may  define  the  vector  la  >  and  i b  >  as 
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This  equation  (22)  can  be  written  as 

\i(V0)l 


'Vs>  * 


hxN(x0’0) 


(23) 


Similarly  using  equation  (E4),  (17),  and  (18),  one  can  write  the 


equation  for  | I  >  as 


I  > 
1  s 


:  Mt 

ey  —I 


eyi(x0,0) 


LeyN(x0*0)J 


(24) 


For  definitions  and  derivations  of  and  whose  columns  are  current 
and  voltage  eigenvectors,  respectively, refer  to  Appendix  (D) .  Note 
that  in  the  above  derivations,  the  energy  stored  in  the  aperture 
junction  was  neglected,  therefore  the  Figures  (8)  and  (9)  represent 
the  zeroth-order  approximation. 

The  first-order  equivalent  circuit  is  now  delivered  for  the  small 
aperture  with  a  MTL  behind  it.  Consider  Figure  (7)  again  with  the 
modification  shown  in  Figure  (10)  with  no  incident  wave  coming  from 
right.  Using  equation  (13)  and  (14)  and  noting  that  the  external 


field  (or  Sc  )  is  zero,  one  can  write 

Tsc+ 


:  »  +i(uu  a  H" 

mx  J  m,xx 


-jwe  a 


„sc+ 


(25a) 


(25b) 


ey  e  y 

Decomposing  electric  and  magnetic  fields  into  their  modal  distributions, 
the  j-th  mode  for  each  field  is 


23 


(31a) 


(31b) 
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e  y  i  yj 

Equations  (30a)  and  (30b)  can  Chen  be  written  in  terms  of  £  and  H  as 


(I  +  j  H  +  j  £)  j  aL> 

(32a) 

(~j  H,  +  j  1)  |aL> 

(32b) 

where  _I  is  an  identity  matrix. 

The  signal  flow  graph  of  equations  (32a)  and  (32b)  is  shown  in 

Figure  (11)  where  ja  >  and  |b  >  represent  the  zeroth-order  equivalent 

s  s 

source  models.  One  may  define  the  transmission  matrix  _T  as 


I*  I+j  H  +  ]E 

and  the  reflection  matrix  R  as 


(32c) 


R  -  -j  H  +  j  E 

Then  equations  (32a)  and  (32b)  can  be  written  as: 


r  n 

—  - 

r  "T 

|b  > 

R  T 

'  ]  a  > 

1  L 

'a.>> 

T  R 

|  0> 

U . 

_  _ 

_ 

(32d) 


(32c) 


Incorporating  the  same  analysis  for  |bR>,  one  obtains  the  complete 
source  free  signal  flow  expressions  as 


f|bT> 


1  a  > 
'  R 


'V 


lVu 


(32d) 


In  the  next  chapter,  this  signal  flow  representation  and  che  aperture 
signal  flow  sources  are  combined  as  the  basic  model  of  aperture 
coupling  to  a  multiconductor  transmission  lines. 


CHAPTER  IV 


BOUNDING  VOLTAGE  AND  CURRENT 
AT  A  TERMINATION  OF  THE 
MULTICONDUCTOR  TRANSMISSION  LINES 


In  this  chapter,  the  procedures  for  calculation  of  maximum  voltage 
and  current  at  a  termination  of  MTL  over  a  small  aperture  excited  by  an 
external  electromagnetic  pulse  is  formulated.  Consider  the  aperture 
representation  to  be  as  shown  in  Figure  (12) . 

Figure  (12)  can  be  represented  by  the  corresponding  signal  flow 
graph  in  Figure  (13).  In  Figure  (13)  T?  and  R  represent  the  trans¬ 
mission  and  reflection  matrices  as  derived  in  equations  (32c)  and  (32d) . 
The  quantities  |ag>  and  |bg>  are  the  source  representations  as  obtained 
in  equations  (17)  and  (18) .  The  quantities  _T^  and  are  termination 
reflection  matrices  as  computed  in  equations  (E15)  and  (E17) .  Finally, 
4^  and  4^  are  the  propagation  matrices  whose  elements  are  exponential 
functions  of  (2)  representing  the  phase  or  time  delays  as 


Diag  [e"jV4,  .  . 
Diag  [e'^lS,  .  . 


»~^SNZ4]  ! 


(33) 


where  g  is  the  phase  constant  for  each  mode.  As  shown  in  Figure  (12) 
the  origin  is  taken  to  be  the  region  of  the  aperture  for  computation 
purposes.  Thus,  [a  (0)>  and  I b_ (0 ) >  represent  the  traveling  waves 
immediately  to  the  right  of  the  aperture  and  |aT(0)>  and  j b ^ ( 0 ) >  denote 
the  traveling  waves  immediately  to  the  left  of  the  aperture. 

One  can  deduce  the  following  equations  by  method  of  signal  flow 
theory: 


26 


Of 
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Figure  (12).  Terminated  MTL  over  an  aperture. 


Figure  (13).  Signal  flow  graph  of  Figure  (12). 
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|aR(0)>  *  T_  |aL(0)>  +  R  |bR(0)>  +  |as> 


|bL(0)>  *  T  |bR(0)>  +  R  [aL(0)>  +  |bs> 


Also  at  terminations  I\  and  r~  one  can  write 

-4  -3 


|aRa4)>  -  4  ]aR(0)> 

I W  *  41  lbR(0)>' 


(34a) 

(34b) 

(35a) 

(35b) 


|bL(-V>  ~h  ibL(0)> 
;ai(-i3)>  ,  ^1|aL(0)>, 


(36a) 

(36b) 


and 


IVV>  ’  I*  I W’ 


(37a) 


|aL(-i3)>  -  r3  ib4(-i3)>  •  (37b) 

Substituting  equations  (35)  and  (36)  in  (37),  one  can  conclude  that 

(38a) 
(38b) 


1  V0)>  m  hhh  |aR(0)> 


|*L(0)>  -  9^  r_3  ^  |bL(0)> 

Having  introduced  equations  (34)  and  (38),  one  can  solve  for  jaR(0)>, 
|b_(0)>,  |aT(0)>,  and  |b  (0)>.  For  this  first-order  formulation, 

A  L  id 

| b^ (0) >  can  be  shown  to  be 

[bL(0)>  -  [(i  -  r3  iL)  -  T(^:  r^1  jT1  -  r) _1  (T  ?l  r3  *  )]  x 

1 IV  -  K41  r;1  !«,>> 


Similarly  |aR(0)>  can  be  obtained  as 


(39a) 


29 


!as(0)>  -  id  -  £  **  h  !*>  -  K41  £  k1  -  «‘1  <1 4  4  4>' 


*  HV  +  K41  Lj1  »[l  -  w 


-1 


lb3>]* 


(39b) 


For  almost  all  practical  problems  of  interest,  one  can  neglect  R 
and  assume  T_  to  be  unity,  meaning  that  there  are  negligible  reflections 
on  the  MTL  due  to  the  aperture.  In  fact,  Kajfez  [12]  has  shown  chat 
the  reflected  amplitudes  are  approximately  one  percent  (or  less)  of 
the  amplitudes  of  the  incident  modes.  Therefore  substituting  R  *  _0 
and  T  *  _I  into  equations  (39a)  and  (39b),  one  obtains  the  following 
equations : 


jbL(0)>  -  II  -  4  r*  4  4  l3  ijrtlv  +  4^4  |.,>] 
|aR(0)>  -  [ i  -  4  r3  ^  %  £4  y1tUJ>  +  4  r3  4  |bs>j 


(40a) 

(40b) 


The  above  two  equations  are  the  zeroth-order  formulation. 

Having  derived  |bT(0)>  and  ]a  (0)>  ,  one  can  use  equations  (36a) 

L  K. 

and  (36b)  to  solve  for  jb^(0)>  and  |a^(0)>  as 

!V0)>  ■  U1  k1 '  k  h  krl['V  +  4  4  4  IV1  (4ll) 

r'1  41  -  4  ^  *Rr1[|b3»  +4^4  i«s>  (4ib) 

Figure  (14)  represents  Che  zeroth-order  signal  flow  graph  for  the 
terminated  MTL  over  an  aperture. 

Having  formulated  simple  expressions  for  |aR(0)>  and  |b^(0)>,  one 
can  calculate  |a(z)>  and  jb(z)>  at  any  point  on  the  line.  Specifically 

for  voltage  at  (.1,)  termination,  one  obtains  the  following: 

*+ 


4 
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|V(z)>  -  My[ !a(z)>  +  |b(z)>] 

|VU4)>  -  Mv[|a(l4)>  +  ibU4)>J 

However,  I  a. C ^4 )  >  and  j b  ( x-4 ) >  can  be  represented  as  follows  from 
Figure  (14): 

|a(*4>>  *  ^  |aR(0)>  (42a) 

|bU4)>  "  I4  1r  |aR(0)>  .  (42b) 

Therefore, 

|V4>  -  |V(i4)>  »  My[I  +  IV]  ^  |aR(0)>  (43a) 

Substituting  (40b)  into  (43a),  one  obtains 

IV  ■  Vi  +  V'4)(i  -  *l  r3  k  4 h  41'1  * 

Mas>  +  I3  ^  |b,>]  (MW 

Using  equation  (  E4  )  and  substituting  from  (43a)  one  can  also  obtain 
any  expression  for  \l,>  as 

•4 

!i4>  -  |IU4>>  -  Mjti  -  rjli  +  I^]"1  M^1  |v4>  (44) 

In  order  to  find  an  upper  bound  for  the  voltage  or  current  at  a 
termination,  |v>  and  jl>  should  be  maximized  in  some  mathematical 
sense.  The  best  mathematical  procedure  for  bounding  vectors  or  ma¬ 
trices  is  the  calculation  of  their  norms.  The  norm  of  a  vector  or  a 
matrix  is  a  single  number  which  gives  both  a  conceptual  and  mathemati- 
can  evaluation  of  the  size  of  the  vector  (or  matrix)  in  the  same  sense 
that  the  modulus  does  for  a  complex  number.  The  norm  of  a  vector  a> 
is  denoted  by  jj  a>  -j  and  satisfies  the  following  relations  [17]: 


II  I*5*  II  >  0 


unless  |  a  >  -  Q_ 


I!  3ja>  j|  *  1 0 1  ||  |a>  ||  for  3  *  complex  scalar 

||  | a  >  +  jb  >  ||  <  ||  |a  >  ||  +  ||  lb  >  || 

||  |  a  >  -  |b  >  ||  >  ||  | a  >  ||  -  ||  |b  >  11 


The  p-norm  of  a  vector  is  defined  as 


a  > 


(-1  ai  I F  +  |a>>r  +  •••  +  la  lP)P,  P  *  1,2, , 


For  the  case  of  p  *  this  norm  becomes 
H  U  >  IL  “  maximum  |a  | 


(45 


(A6 


Similarly,  the  associated  norm  of  a  matrix  A  is  denoted  by  ||  A  || 
and  satisfies  the  following  relations: 

j|  A  ||  >  0  unless  A  »  0 

||  3  A  i)  *  |  S  |  I)  A  ||  for  3  3  complex  scalar 

II  A+ B_||  ^  II  A  ||  +  ||  _B  j|  (triangle  inequality) 

II  A  B  |!  <  |f  A  ||  ||  B  || 


The  norm  of  a  matrix  for  special  values  of  (p)  is  defined  as 

II  *  II  _  E  I-  I 


U, 
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vector  is  greater  or  equal  to  its  corresponding  infinity-norm  since 


Dower  : 


<  b  |  b  >  <  <  a|a>  for  physically  realizable  systems 


or 


<  b|b>  »  <  a|r_+_T|a>  <  <  a|a> 


C51) 


Similarly,  for  any  eigenvector  |u^>  of  matrix  with  eigenvalue 
X ^  one  obtains 

<  ui|r+r|ui >  -  <  ui|ui > 

But  according  to  equation  (51) 

<  ui|r+r|ui  >  <  <  uiju1  >  . 


Therefore 

(Xf)  <  ui  |  ui  >  ^  <uilui> 

or 

(Xt<  1) 

Particularly 


X  <  1 
max 


or 


v  Xm  _  „  ^  1 

max 


(52) 


The  above  equation  is  the  definition  of  the  tvo-norm  of  matrix  and 
thus 

II  I  !|  <  1  (53) 


At  this  3tage,  one  can  start  analyzing  equation  (43b)  by  taking  the 
infinity-norm  of  |V,  >  as 


II  |v4>  II  -  Hu* 


las>+WLlb3>^l 


Using  the  triangle  inequality  for  norms  and  other  relations  as  intro¬ 
duced  previously,  the  above  equation  is  easily  shown  to  be 


II  h4»ll.<lla»ll.tll i II. +114 ll„ -’ll  1*11  Jit 


-1 


III  Is >H  +  ll*,H  llrjl  114,11  II  lb. 

*3  00  L»  00  -^ao  Loo  - 


>11  ] 


Invoking  equations  (48)  and  (43b),  the  above  inequality  can  be 
expressed  as 

-•n  11  1 


II  |V4  >  ||  <  IlH^II  til  I  II  +  III 


,  •  11J4  11,  J  11  inn, 

2  2  2  1 .  U  r  <6  A  r  $ 

1  1 1  -  L“3— -L~R~4~R  2 


HI  IS  "I  L*  Hi  UJllJIJliJIJI  Ib,>||  1  (54) 


'2‘  "^L  V'~3  "2  ”~L2 


Using  equation  (33) ,  one  obtains 


(55) 


Substituting  equations  (55)  and  (53)  into  equation  (54),  one  has 

(  ||  I||2  -  1  ) 


Is 


1  - 


4  T  <4  ST, 4_ 
-  L-3-  L-St-4-R 


las> 


II  Ibg> 


iU 


The  norm  in  the  denominator  may  be  expanded  as 


4  r„* 


(56) 
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◦r  by  substitution 

II  !LI3WA»,  < 1  • 


But  in  order  to  make  the  righthand  side  of  equation  (.56)  a  finite 
quantity,  one  has  to  mathematically  introduce  the  practical  losses 
associated  with  each  mode  travelling  on  the  MIL  as 


' 

4 

•  a 

Diag 

-vl2  -yNz 

e  ,..,,e 

$ 

~L 

Diag 

+Ylz  +yNz 

e  ,...,d 

4 

(57) 


where 


Y±  *  «i  +  j  Si 


(58) 


is  a  decay  constant  for  each  mode. 


Therefore  by  (,47b) 

II  $ 


max  .  "ai^4. 
i  Ce  5 


(59) 


L  2 


max  ,  ”aii3. 
i  (e  5 


where  (l^)  and  (£3)  are  right  and  left  termination  distances  from  the 
aperture  as  depicted  in  Figure  (12).  The  above  equation  can  also  be 
written  is 

"  (aiaiin  )J,4 

(60a) 

-Knin^ 


e 


II  A  II  _  - 


Equation  ( 56 )  may  now  be  expressed  as 


-  e'(a^ 


dL) 


'V 


b  > 

1  3 


where 


°R  '  +2aimin14 

%  m  +  2aW3 


If  is  defined  as  the  sum  (cR  +  0^) ,  then 

*  2oimin^4  +  "  2ainiiniT  (60b) 

thus 


|v4>  || 


H,l 


Cl  -  S*) 


III 


is>l!2 


ib«  >iU 


(61) 


The  quantity  oT  may  also  be  interpreted  to  include  losses  in  the  termi¬ 
nation  corresponding  to  non-unity  bounds  on  the  |£^>  and  to  include 
radiation  loss  both  back  through  the  aperture  and  in  coupling  to  higher 
order  modes. 

It  is  computationally  advantageous  to  replace  ||  H,,  with  the  two-norm 

which  is  consistent  with 


!!  IV  II  *  II  IV  H, 

00  4 


(62) 


This  may  be  seen  by  determining  the  two-norm  which  by  equation  (D27 ) 
and  (D31)  is 


(63) 


-1 


where  the  X/a  are  the  eigenvalues  or  the  MTL  matrix  L  which  is  the 


inverse  of  induction  matrix  of  the  MTL.  Also  using  the  triangle  in¬ 


equality,  one  obtains 


as>  H2*  [Js  Hcmxfhx>  !l2+  4  II  “cey|e7>  11 2  J 


and 


II  lbs>  Hj*  tJs  H“cmxlbx>  H2+  **  H'ceyley>  U2J 


Since 

||  c  |h  >  ||  -  ||- c  |h  >  || 

11  mx1  x  11 2  mx 1  x  2 


j |  —  c  |  e  >  ||  ■  ||  c  | e  >  || 

11  ey  y  2  1  ey  y  2 


equation  (61)  may  be  written  as 


2 

a 


mx 


c  |  e  >||  ] 

ey 1  y  "  2  J 

(64) 


Finally,  since  the  vector  |V>  defines  the  voltage  on  a  MTL  with  respect 
to  ground,  the  maximum  voltage  of  interest,  in  other  words  the  voltage 
between  any  two  lines,  is  bound  by 


AV  <  2  |  V,  > 

1  max  11  1  4 


Using  the  above  inequality,  equation  (64)  may  be  written  as 
4 


AV  < 


“*  Cl  -  e‘aT) 


/max (  )  [  ||  c  |h  >  ||  +  l|  c  |e  >  i|  ] 

v  i  "  mx1  x  11  ?  1  ey1  y  "2 


(65) 


which  is  the  desired  upperbound. 

Substituting  for  c  and  c  from  equations  (13)  and  (14),  equation 

w  fflY  A  V 
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(65)  can  be  expressed  as: 


Uv 


'max 


4u 


_  e-<JT 


+  £ 


let!  j  E  3C* 

e  bound  y 


-  ^ 
1  y 


(66) 


Having  determined  a  maximum  bound  on  the  termination  voltage,  one  can 
calculate  a  maximum  bound  on  the  termination  current  using  equation 
(43a)  and  (43b)  to  simply  obtain 


K>M  I!  !V!I2 


|MVI 


(67) 


or 


AV  Mt 

1  'max  "  —I  n2 


max' 


—V  I!  -i 


(68) 


Invoking  equations  (D28)  and  (D32),  one  can  calculate  ||m  Jlj  in  the 
same  way  as  ||  My  Co  obtain 


M^||^  »  J max  (eigenvalue  of  H^.1  M^)  *  J max(  ) 


(69) 


Therefore  equation  (68)  can  be  written  as: 
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or 
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vi  Vj 

maxC —  )m±nC  -r^ 

\i  Ai 


C71) 


Therefore  by  equations  (66)  and  (71)  one  has  a  bound  on  both  voltage 
and  current  at  a  passive  termination  located  at  some  distance  from  the 
aperture.  Although  the  above  formulations  all  have  been  based  on  the 
fact  that  the  termination  is  to  the  right  of  the  aperture  at  a  distant 
(£4),  equations  (66)  and  (71)  provide  a  general  frequency  domain  upper 
bound  at  any  point  on  the  MTL.  It  should  be  noted  that  the  losses  due 
to  the  aperture  or  termination  can  be  taken  into  account  by  modifying 
oT  as  previously  suggested. 

The  next  chapter  transforms  these  results  to  transient  domain,  and 
includes  the  important  multiple  reflection  phenomenon. 


CHAPTER  V 


TRANSIENT  ANALYSIS 

The  previous  discussions  were  based  entirely  on  the  frequency 
domain  analysis  where  each  mode  was  treated  as  steady  state  sinusoidal 
function  of  time.  In  this  section  the  transient  (time-domain)  analysis 
is  considered  together  with  the  multiple  reflection  phenomenon. 

A  quasi- TEM  analysis  of  the  MTL  is  discussed  in  Appendix  (E) .  It  is 
assumed  that  this  quasi- TEM  waveform  is  non-dispersive,  in  other  words 
an  arbitrary  waveform  is  transmitted  by  each  mode  without  distortion, 
and  the  waveform  at  the  distance  (2^)  is  a  replica  of  the  transmitted 

A4 

waveform  delayed  in  time  by  (  —  )  where  V£  is  the  velocity  of  the  mode. 
Mathematically,  this  is  expressed  as 

a^^.t)  -  a^Co.t  -  —  )  (72) 


The  signal  flow  graph  for  the  time  domain  is  given  in  Figure  (15) 
where  terminations  are  located  at  distances  (2^)  and  (2^)  to  the  right 
and  left  of  the  aperture.  If  the  ^  matrix  of  the  previous  chapter  is 
given  by 


h 


.-38*4 


(73) 


then 

|aR(24,t)  >  -  ja^o.t  -  t4)  > 


(74) 


where 


(75) 


The  above  formulation  implies  that  all  the  modes  have  the  same  velocity 
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and  consequently  the  sane  phase  constant.  Therefore  it  takes  units 
of  time  for  all  modes  to  impinge  on  the  termination.  On  the  other  hand 
if  the  matrix  is 

K  -  [5.,  e"^1*4]  (76) 

— R  ij 


tl >  -  |  [a^Co.t  -  ri)]  > 


where 


1  vi 


In  this  case,  each  mode  has  its  own  characteristic  time  constant  (t^), 
and  the  formulation  is  much  more  complicated.  In  the  following 
analysis,  there  are  three  major  assumptions  that  should  be  kept  in 
mind:  1)  The  propagation  matrices  and  are  taken  to  be  constant. 
Mathematically,  it  is  written  as 


U! 


-je ik  ) 


where  3  is  defined  as 


«  •  7 


Also  che  variable  (t)  corresponding  to  the  transit  time  from  (-i^)  to 
(+4^)  is  defined  as 


*3  +  l4 
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2)  Terminations  are  assumed  to  be  passive  and  resistive  which  implies 
that  the  elements  of  termination  matrices  £^  and  £^  are  real  and 
constant;  3)  The  medium  permittivity  ( e )  is  taken  to  be  that  of  free 
space  (eQ)  so  that  all  the  modes  travel  with  speed  of  light  (v»c). 

The  above  three  assumptions  imply  the  worst  case  conditions  for  the 
whole  system  in  terms  of  the  voltage  built  up  at  terminations  as  a 
consequence  of  having  the  highest  possible  speed  for  each  mode 
traveling  on  the  MTL.  This  directly  relates  to  minimal  mode  attenua¬ 
tion  between  successive  multiple  reflections.  Also  as  implied  in 
Figure  (15),  the  aperture  transmission  and  reflection  matrices  are 
assumed  to  be  unity  and  zero  as  given  by 


1*1 
R  -  0 


(80) 


This  assumption  simplifies  the  formulation  and  calculation  of  the  pro¬ 
blem  and  due  to  the  very  small  effect  of  £  and  R  for  problems  of 
interest  [12J,  it  is  a  good  approximation  to  the  exact  solution. 

In  Chapter  IV,  the  termination  voltage  was  calculated  as 

IV  -MyU  +  V  4  |  aR(0)>  . 

In  the  time  domain,  4^|aR(0)>  is  replaced  by  a  series  of  incident  and 
reflected  waves  constituting  the  multiple  reflection  phenomenon  as 

!V  >*M±  +  l41  (  lasCt“V>  +  —3^4 1  a3 (t  “  ”4  ~  “T }  > 


as(t- 


-n(27))>  + 


4- 
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!3ib3Ct  +  T4~2T)>+  2AI3^3Ct4.T4-2C2T))> 

'l 

+  ...  +  (r3r4)ar3jbsct  +  T4-nC2T))>  + ...  j  (8i) 

The  above  equation  can  be  best  understood  by  observing  Figure  (16) . 

At  time  (t),|as(t}>  travels  (T^)  units  of  time  to  reach  the  right 
termination.  After  being  reflected  by  r^,  the  wave  travels  (t)  units 
of  time  to  reach  the  left  termination,  and  again  is  reflected  by  F_3 . 
Finally,  it  travels  Or)  units  of  time  to  impinge  on  the  right  termina¬ 
tion.  This  process  continues  indefinitely  as 

CI^)11  i  as(t-T4-nC2r))>. 

Similarly,  one  can  deduce  that  the  process  for  | bg (t ) >  takes  the 

form 

CI^)11^  !  bg0t  +  T4-n(.2T))>  . 

Equation  (81)  can  now  be  simplified  to 

00 

|V4>- r^]  [  l  (r3r4)n!as(t-T4-n(2x))>  + 
n-0 

00 

l  (I3£4)nl3|bs(t  +  T4-2(n+l)T)>]  (82) 

n-0 

Due  to  the  losses  of  the  aperture,  the  terminations,  and  the  MTL,  the 
above  equation  may  be  modified  to  include  such  losses  as 

00 

I'V-SMI  +  IJI  l  e"0n(I3^)I1|as(t-T4-a(2T))>  + 

n«0 
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I  e"PnQUI 4)nI3ibaCt  +  T4-2Cn+l)T)>J  (83) 

n»0 

where  (e  pn)  takes  into  account  the  system  loss  for  each  time  that  a 
wave  travels  twice  the  length  of  the  MTL  and  passes  over  the  aperture. 
In  order  to  simplify  the  computation  of  equation  (83) ,  one  can  use  a 
pole  expansion  method  of  the  Singularity  Expansion  Method  (SEM)  [18]. 
This  expansion  may  he  written  for  simple  poles  of  a  general  function 
F(t)  as 

m 

F(t)  -  If,  eSJC  u  ( t )  (84) 

j-1  3 


where  the  s^'s  are  simple  poles  of  the  transform  of  F  (t) ,  and  f^’s  are 
residues  due  to  each  pole.  Using  the  format  of  equation  (84)  in 
equation  (83),  one  can  write 


(,t-T4-n(.2T))>  -  I  |aSj(0)  >esj(t“  T4_n(2T))u(t-T4-n(2T)) 


j-1 


(35) 


m 


|bs(t  +  t4-2(n+l)T)>=  i  j  bs,  (0)  >  e  3  v  4 

j-1  2 


s-j  (t  +  t,  -  2(n+l)t) 


u(t  +  -  2 (n+l)t) 

(86) 


where  (s^)  Is  Che  propagation  constant  for  each  external  made  in  the 
function  expansion.  If  equations  (85)  and  (86)  are  substituted  in 


equation  (83),  one  obtains 

|V,>«  I  eS:(C"T4-n(2T))e"0n  + 

•  '  '  1  n=0 
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I  |bs.(0)>  + 


4-1  3  n-0 

where,  due  to  Che  unit  step  functions, 

t  “  t. 


(87) 


N  *  Integer  [ 


2t 


] 


t  +  x,  -  2t 

M  -  Integer  [  - — -  ] 


(88a) 


(88b) 


for  N  and  M  non-negative.  If  N  or  M  is  negative,  the  corresponding 
series  is  deleted. 

In  order  to  obtain  an  upper  bound  on  the  voltage  at  the  termination, 
one  should  take  the  two-norm  of  equation  (87).  Using  the  matrix  and 
vector  norm  properties  introduced  previously  along  with  the  triangle 
inequality,  one  obtains 


1  ; 
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|v4> 

2 

2 

1  + 

14 
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j-i  L 

lasj  (0)> 
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l  ^ 

n*0 

M 


s-j  (t  -  T4  -  n(2x))g-on 


bSj  (0)> 


y  esj(t +  -  2(n+l)x)e-pn| 

'n-0  ‘ 


(89) 


If  Sj ,  the  propagation  constant,  is  decomposed  into  its  real  and 
imaginary  part  as 


s(j)  *  "n(j)+j  5U) 


tei 


f  e3j(-t“T4_tl(2T))e'pn 
n”Q 


<  ?  e"nj(t'T4*Q(2T))e~pI1 

n»0 


and 


l  eSjCt  +  x4  -  2(n+l)T)e-pn 
n*0 


^  y  e-njCt  +  T4-2Cn+l)T)e-pn 
n*Q 


The  quantities  (-r^)  and  (£^)  represent  the  decay  and  phase  constants 
respectively  for  each  mode. 

One  can  use  the  finite  geometrical  series  to  simplify  the  above 
expressions  as 


y  -nj(t-  T4-nC2f))e-pn  ,  e-nj(t-t4-2NT)e-pN. 
n**0 

”  1  _  e_(N+1)  (.2*1  jx  -  p )  ” 

_l_e-(2T1jT‘p) 

and 

y  e-nj(t  +  T4-2(n+l)T)e-pn  a  e-nj(t  +  T4-2T-2Mx)e-pM 

n*0  _ 

_  e-(MH)  (2njX  -  p)  ~j 

l-e'C2njT  "  0) 


(90) 


(9i: 


Equations  (90)  and  (91)  hold  only  when 


2n,x  >  p 
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Ac  this  stage,  one  can  start  computing  a  numerical  upper  bound  on  the 
voltage  at  the  termination. 

The  maximum  value  for  both  e  T4  an(j  e  +  T4  2t  -  _Mr)  is 

unity  which  occur  approximately  when  (N)  and  (M)  are  replaced  by  their 

respective  values  given  in  equations  (83). 

In  order  to  find  an  upper  bound  for  the  remaining  terms  in  equations 
(90)  and  (91) ,  one  can  set  M  equal  to  N  and  compute  the  values  of  N  for 
which  the  remaining  terms  are  maximum.  That  is 


-ON  li-.-Wg'l'-’ 

8  U-e-<2rl3T-a> 


After  caking  the  derivative  and  solving  for  N  denoted  by  N' ,  one  ob¬ 


tains 


£nC2n:)T/o) 

S'  *  (2rijT-p) 


The  above  value  for  M'  maximizes  the  equation  which  can  be  verified  by 
showing  that  the  second  derivative  is  negative.  Since  S'  must  be 
greater  chan  or  equal  to  zero  for  the  correct  solution,  one  must  com¬ 
pute  (95)  and  replace  N''  with 
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N  -  max  [0,Nr] 
max  J 


(96) 


Substituting  N  in  aquations  (9Q)  and  (91) ,  one  obtains 
niflx 


l  e"nj  (-c  "  T4  "  nC2T))e-pn  <  • 
n*0 
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f  2r.4t  ~ 
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(97a) 


and 
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-pn  ^  • 
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2rux 

An  (—■*-) 

0  +  2p,T 
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» 

—  — 

- 

(1  --=-*= — > 
2n^x 


1  -  e 


-(2n.x-p)  f 

J  j 


(97b) 


for  N'  greater  than  zero  and  unity  otherwise.  For  the  case  when 
(2hjT  *  p),  equations  (92)  and  (93)  are  used  and  the  results  can  be 
shown  to  be 


N 


I  e'n3 

n-Q 


-rn(t  -  t4  -n(2x))e-pn  <  \ 


(i).*"15.  P<1 
1  ,  P>1 


(98) 


The  above  results  also  hold  for  the  summation  over  M. 
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Equations  (,97a)  and  (,97b)  may  be  simplified  notationally  by  defining  a 
new  variable  y  as 

y  -  2n^t  . 

Hence,  the  above-mentioned  equations-  can  be  written  as 

?  e-njCt-T4-nC2T))e-on  K 
n*0 


pln( P/  )  " 

- 

exp 

p  + - 2- 

Y  -  P 

»  — 
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1 


(p/Y> 

e(p-Y) 


(99) 


and  similar  simplified  result  holds  for  summation  over  M.  If  (p)  is 
zero,  the  system  has  ideally  no  internal  loss  and  the  upper  bound  of 
the  voltage  at  the  termination  will  be  maximum.  Thus  the  worst  case 
condition  suggests  no  loss.  Substituting  zero  for  (p)  in  equation  (99), 
one  obtains 


f  e'nj  (t  "  T4  -n(2t))e-on 
n*0 

f  e_n j +  T4  “  2 (n+1) x) e-pn 
n«0 


1 

(1  -  e”Y) 


(100) 


Finally,  substituting  equation  (.100)  in  equation  (89),  one  obtains  the 


voltage  upper  bound  as 
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(101) 
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With  Che  simplif ications  of  equations  (.53)  and  (63),  equation  (101) 
becomes 

I  r 

'  "  !*„»)>  II,  +  II  |b,,(0)> 
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U  'J 


9 J  '  2j 


-2n,  t 
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(102) 


As  discussed  in  chapter  IV  ,  the  maximum  voltage  between  any  two  lines 
is 
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AV  |  2 

|vA> 

|  max 

1  4 

Therefore 


A  V 


max 


*/T^ 
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l»Sj  (0»  II  2+  II  I  bSj  (0)>  ||, 
___ 

(1  -  e  2  ) 


(103) 


One  can  use  equation  (70)  to  obtain  an  upper  bound  on  the  current  as 


AV 


'max 


max 


/ 

v  max,v  , 


,  /  v  \  min  /  v  , 
i  i 


(104) 


If  |ag^(0)>  and  [bSj(0)>  are  replaced  by  their  equivalent  express¬ 


ions  in  (hx>  and  |ey>,  then  equation  (103)  becomes 

<  4  /“f  (^)  {  l 
max  v  i  i  L 


cmx(0)  hx>,!2  ''  cey  (0)  |ey  > 


j-l 


-2n.t 
(1  -  e  2  ) 


(105) 


In  time  domain,  c  and  c  can  be  expressed  as 
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Cmx--uVxs£  [HxSC’(t)J 


and 


'ey  *  e*e£ 


(106a) 


(106b) 


where 

H  SC“(t)  -  (h/c~)  1  f,  eSJC 
j»l  J 

E/C'(t)  -  CEyac->  ?  t.  e3J  (t) 
7  y  j«l  3 


(107a) 


(107b) 


Taking  derivatives  from  equations  (107a)  and  (107b),  and  substituting 
in  equations  (106a)  and  (106b)  respectively,  equation  (105)  can  be 
written  as 


i 
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m,xx 


bound 


ae  bound 


(108) 


for  a  lossless  MTL. 

Equations  (108)  and  (104)  suggest  a  computational  method  for  upper 
bound  signal  levels  in  time  domain  at  a  termination  of  MTL  behind  an 
aperture-perforated  conducting  screen. 

In  the  next  chapter,  the  resultant  equations  will  be  used  to  com¬ 


pute  an  upper  bound  signal  level  for  a  specific  problem  both  in 
frequency  and  time  domains. 


CHAPTER  VI 


COMPUTATIONS  OF  SIGNALS  UPPER  BOUNDS  FOR  A 
PARALLEL-PLATE  TWO-CONDUCTOR  TRANSMISSION  LINE 


In  Chis  chapter,  the  results  of  previous  discussions  along  with 
Appendices  are  used  to  compute  an  upper  bound  for  the  voltage  and  the 
current  at  a  termination  of  a  parallel-plate  two-conductor  transmission 
line  located  behind  a  circular  aperture.  Hopefully,  this  example  will 
help  the  reader  understand  fully  the  concepts  and  the  procedures  used  in 
previous  discussions  in  computations  of  upper  bounds. 

The  four  basic  equations  which  fully  describe  upper  bounds  on  the 
voltage  and  the  current  in  both  time  and  frequency  domains,  as  derived 
in  chapters  IV  and  V,  are 
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C104) 


The  parameters  used  in  this  particular  problem  are  the  same  as  the  ones 
introduced  in  an  example  by  Kajfez  [12]  which  facilitate  comparison  of 
results. 

There  are  a  few  concepts  that  should  be  mentioned  before  proceeding 
to  actual  computations.  By  transmission  line  theory,  the  induction 
matrix  L  is  related  to  the  induction  coefficient  matrix  K'  as 


where  c  is  the  speed  of  light  in  vacuum.  The  matrix  K''  describes  the 

capacitance  relationships  of  the  homogeneous  MTL  and  is  evaluated  for 

the  system  filled  with  a  vacuum.  This  matrix  depends  only  on  the 

geometry  of  the  system.  The  quantities  X^'s  appearing  in  all  the  above- 

mentioned  formulas  are  the  eigenvalues  of  the  L  ^  matrix,  and  the  v^'s 

2 

are  the  corresponding  modal  velocities  related  to  the  eigenvalues  (-3^  ) 
2 

of  the  matrix  (.-»  LK)  as 


u 


where  K  is  the  same  matrix  as  K'  except  that  it  is  evaluated  with  the 
true  values  of  permittivity  of  the  medium  contained  in  the  system. 

The  parallel-plate  two-conductor  transmission  line  is  located 
behind  a  circular  aperture  of  diameter  d  ■  2  cm  which  is  positioned  at 
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z  ■  0  as  shown  in  Figure  (17).  The  transmission  line  is  terminated 
by  matrices  and  representing  passive  terminations  located 
distances  of  5m  and  7m  from  the  aperture  respectively.  The  cross 
section  of  the  transmission  line  is  shown  in  Figure  (.18).  There  are 
three  layers  of  dielectrics  separating  two  strip  conductors  of  width  W. 
The  thickness  of  dielectrics  are  denoted  by  h^,  f^,  h^,  and  are  assumed 
to  be  much  smaller  than  the  width  W. 

The  computations  of  both  the  sources  and  the  induction  matrix  require 
knowledge  of  the  quasi-static  modes  on  the  MTL.  In  order  to  obtain  the 
electrostatic  field,  a  unity  potential  is  applied  to  conductor  a,  while 
conductor  b  and  the  shield  are  held  at  zero  potential.  Then,  the 
normalized  fields  in  three  regions  are 


Similarly,  a  unity  potential  may  be  applied  to  conductor  b  while  conduc¬ 
tor  a  and  the  shield  are  held  at  zero  potential.  Then,  the  normalized 
fields  in  three  regions  are 


(110) 


Computation  of  the  elements  of  the  induction  coefficient  matrix  K  is 
accomplished  by  integrating  the  electric  flux  through  the  closed  surface 
Sa  around  the  conductor  a  as 
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Substituting  these  parameters  in  equations  (112),  C113),  and  (114),  one 
obtains 
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where  (£^)  is  the  permittivity  of  vacuum.  If  all  three  dielectrics  are 
assumed  to  be  air  (  »  e0  **  £^  ■  £q) ,  then  the  induction  coefficient 

matrix,  denoted  by  Kf ,  becomes 


15 

-10 

"  e0 

-10 

15 

Thus  L  ^  can  be  obtained  as 
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where  is  the  permeability  of  vacuum.  The  corresponding  eigenvalues 
of  l  1  can  be  computed  as 


“O 


2  2 

The  eigenvalues  (_-  8^  )  of  Che  matrix  B_  ■  (.-w  LK)  and  the  correspond¬ 
ing  eigenvectors  can  be  obtained  as 


cf)2 


1 

n 


and  the  corresponding  modal  velocities  are  computed  using  the  relation 


where  c  is  the  speed  of  light  in  vacuum. 


To  determine  the  sources,  the  previous  results  must  be  extended  to 
compute  the  aperture  modal  function  !e^>  anc*  |hx>. 

3y  equation  (D31),  the  matrix  tty  which  transforms  a  power  wave  to  its 
corresponding  voltage  value  is  obtained  by  computing  the  voltage  eigen¬ 
vectors  |<(k>  given  by  equation  (D27).  Substituting  the  particular 
values  of  ,  v^,  and  in  equations  (D27)  and  (D31)  one  obtains 
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Similarly,  by  equations  (D29)  and  (D32),  one  can  compute  the  matrix 
Mt  which  transforms  a  power  wave  to  its  corresponding  current  value  as 
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By  definition,  the  electric  field  of  the  nth  mode  traveling  in 
positive  (zl  direction  is 


iTCx,y,z)  -  an  V  e"jSnz  (115) 

where  is  normalized  for  unit  power. 

By  selecting  aQ  =  1,  one  has 

if (x,y,0)  *  e£  (x,y> 


The  corresponding  voltage  vector  for  (aQ  -  1)  and  (bQ  *  0)  is  obtained 
from  equation  (E3)  as 


V  >  -  j$  >  = 

1  n  1  vn 


an 

b 

bn 


If  the  potentials  on  the  two  conductors  are  selected  equal  to  6  and 

an 

4>kn  as  shown  in  Figure  (19),  then  the  modal  function  en  of  the  nth  mode 
can  be  obtained  as 


ev  »  i 
yn  an 


(-E2A}y  +  ^bn  (-E2bK 


"bn 


(116) 


Extending  this  result  to  both  modes  (n  ■  1,2),  one  obtains 


!*v> 


'  0. 


hl  ! 

b2  i 


which  gives 


e  > 
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-307.0  " 
+118.5 


for  the  problem  of  interest. 

The  modal  function  hQ  ( _x,y )  is  equal  to  the  magnetic  field  H^Cx.y.O) 
inside  the  transmission  line  when  the  current  on  conductor  (a)  and  (b) 
are  respectively  selected  to  be  ^an  and  4i^r  as  shown  in  Figure  (20) . 

The  currents  are  assumed  to  be  uniformly  distributed  over  the  con¬ 
ductor  surfaces.  An  elementary  computation  using  Ampere's  law  and 
magnetic  flux  conservation  gives  the  following: 


xn 


♦an  hl  *  ♦bn<hl  +  h3) 
W(h1  +  h^  +  h3) 


(117) 


where  y  and  tp,  are  obtained  from  the  matrix  which  transforms  a 
an  bn  — j. 

power  wave  to  its  corresponding  current  value.  The  vector  |hx>  for 
the  problem  of  interest  becomes 

"0.814  “ 

|  h  >  *> 

X  -0.421 


The  external  short  circuit  fields  are  needed  to  complete  the 
bounding  process.  For  an  incident  wave  with  transverse  magnetic  (TM) 
polarization  as  shown  in  Figure  (21) ,  one  has 


EySC~  -  2Eq  sinO,  HXSC_ 


-2E, 


sina 


(113) 


The  incident  parameters  chosen  by  Kajfez  [12]  are 


Figure  (20) .  Computation  of  the  modal  magnetic  field  h 


Figure  (21). 


Transverse  Magnetic  (TM)  polarization 
incident  on  the  aperture  [12]. 
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0  -  45°  ,  a  -  30°  ,  Eq  -  100  KV/m>  nQ  ■ 

In  order  to  obtain  a  numerical  value  for  equation  (66)  one  must 
compute  the  norms  and  absolute  values  for  the  following  quantities: 
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Assuming  (ai)min  *  '100ft  ’  one  obtains 
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-CaiW100ft) 
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or 


Vatin"  3-777  *  10 


-3 


and  by  equation  (60b),  becomes 
JT  -  90.652  x  10~3 


Now,  one  can  substitute  all  the  available  data  into  equations  (66)  and 


(71)  to  obtain 
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IavI  <  (ui)  (191.4  x  lo-9)  V 
1  1  max 

I  I  |  <  (w)  (3.287  x  10~9)  A 

If  the  dielectric  is  replaced  by  a  vacuum  and  the  incident  angle 
is  set  (0  =  90°,  a  =  90°)  for  maximum  short  circuit  fields,  the  upper 
bounds  are 

| AV I  <  (oj) (4.76  x  10“7)  V 
1  1  max 

and 

III  <  (ui)  (7.06  x  10~9)  A 

This  forms  a  useful  frequency-domain  bound  for  the  problem  presented  with 
a  voltage  less  than  10  Volts  for  frequencies  below  3.3  MHz. 

For  comparison,  this  problem  has  been  solved  exactly  for  open- 
circuit  terminations  on  the  MTL.  For  such  a  case 


^3  =  h  "l* 

Determining  the  |a  >  and  |b  >  of  equations  (43b)  from  (13),  (14),  (17) 
and  (18) ,  the  results  were  computed  and  are  plotted  in  Fig.  22  along 
with  the  bound.  A  modified  bound  is  also  plotted  which  represents  the 
actual  bound  of  the  particular  problem.  The  difference  in  bounds  is 
4.1  which  seems  slightly  unreasonable  until  the  bounding  approach  is 
examined.  A  factor  of  two  arises  in  the  bound  to  account  for  a  dif¬ 
ferential  mode  which  does  not  occur  in  the  case  considered.  The  tri¬ 
angle  inequality  used  in  the  bound  of  | ag>  and  ]bg>  contributes  another 
1.5.  The  product  of  the  2-norms  of  M^  and  | e^>  versus  the  <»-norm  of 
(My|ey>)  contributes  a  1.29  factor.  A  small  contribution  also  occurs  due 
to  some  of  the  neglected  loss  terms.  In  light  of  these  observations,  the 


resultant  bound  is  very  reasonable. 
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For  Che  Transient  Analysis,  the  external  incident  fields  given  by 


equation  (118)  should  be  modified  in  order  to  take  into  account  the 
time  variations  as 


E  SC_(t)  =  2E  F(t)sin0,  H  SC,(t)  *  - -  F(t)sino  (119) 

y  u  x  n  Q 


where  for  EMP  a  suitable  function  is 


F(t) 
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Kaj fez's  date  [12]  are 
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Equation  (120)  is  in  the  format  of  equation  (84)  with  |  f  |  =  l(j=l,2), 
=  a,  and  *  6«  Using  equations  (108)  and  (104)  for  the  vacuum 
filled  lines  with  Kaifez's  data,  one  obtains 


[av|  <  (1.8959)  V 

'max 

|l|  <  (0.0281)  A 

'max 

where  |H  SC_|  and  |E  SC  |  in  equation  (108)  correspond  to  the  ex- 
x  y 

pressions  of  equation  (119)  without  the  F(t). 


Kajfez  [12]  has  obtained  a  result  for  this  problem  with  moderately 


mismatched  terminations.  He  has  obtained  a  maximum  voltage  peak  of 
about  (0.27)  volts  for  conductor  (a)  and  (0.25)  volts  for  conductor  (b). 
Although  at  first  glance,  it  may  appear  that  the  result  obtained  in 
this  discussion  (1.9  volts)  is  not  reasonable  at  almost  ten  times  that 
of  Kajfez,  it  should  be  noted  that  for  upperbound  computations  it  is  in 
reasonable  agreement.  The  reason  being  that  Kajfez  has  not  computed 
worst  case  conditions  and  thus  in  solving  the  problem  he  has  not  de¬ 
viated  from  the  actual  parameters  regarding  the  termination.  With  his 
moderate  mismatch,  multiple  reflection  phenomenon  have  not  affected 
the  bound  due  to  the  loss.  Also,  in  computation  of  various  variables, 
he  has  used  the  true  permittivity  of  the  system.  In  this  discussion, 
the  use  of  multiple  reflection  phenomenon  has  had  a  negligible  effect  on 
the  bound.  At  several  stages  of  computation,  other  matrix  norm  in¬ 
equalities  such  as 

II  AB  |  <  |  A  «  ||B  || 

and 

II  A+  B||  <  I  A  |  +  |  B  I 

have  been  used  extensively.  These  inequalities  have  contributed  twice 
or  more  to  the  exact  values,  and  therefore  | AV J  may  be  reasonably  as 


1 AV|  <  0.95  volts 
1  'max 

In  the  transient  analysis,  the  velocity  of  all  the  modes  have  been 
assumed  to  be  the  velocity  of  light  in  a  vacuum  which  is  actually  larger 
than  the  true  velocities  of  the  modes  in  their  respective  media.  In 
addition,  the  upperbound  has  had  a  factor  of  two  to  account  for  the 


possible  potential  difference  between  wires  and  not  just  to  ground. 
Taking  all  these  factors  into  account,  one  can  readily  observe  that  the 
(1.9)  volts  for  a  bound  on  the  voltage  is  indeed  in  reasonable  agree¬ 
ment  with  Kaj fez's  calculations. 


CHAPTER  VII 


I 


UPPER  SIGNAL  BOUNDS  FOR 
ADDITIONAL  CANONICAL  PROBLEMS 

A  more  complete  understanding  of  the  bounding  procedures  will  be 
presented  in  this  chapter  in  connection  with  the  upper  bounds  of  three 
additional  canonical  problems.  These  problems  are  two  parallel  thin 
wires  behind  an  aperture,  parallel  thin  and  thick  wires  behind  an 
aperture,  and  a  wire  between  an  aperture-perforated  parallel-plate 
transmission  line.  The  latter  two  represent  a  wire  at  the  surface  of 
a  cable  bundle  and  a  wire  between  bulkheads  respectively. 

As  in  the  previous  chapter,  we  use  the  bounds  of  equations  (66), 
(71),  (108),  or  (104)  as  appropriate.  Let  us  first  concentrate  on  the 
voltage  represented  by  equation  (66)  below: 


+  € 


bound 


I  lie  >||  ] 

y  2 


(66) 


Equation  (66)  was  actually  obtained  from  an  upper  bound  on  the  power 

along  the  transmission  line.  The  square  root  is  the  upper  bound  on  the 

transformation  from  power  waves  to  voltage  waves.  The  factor  of  four 

accounts  for  the  wire  to  wire  voltage  rather  than  just  the  wire  to 

ground  voltage  in  addition  to  the  total  termination  voltage,  which  can 

-CJy 

not  exceed  twice  the  incident  voltage  wave.  The  quantity  (1  -  e  ) 
accounts  for  the  total  multiple  reflections  with  representing  the 
loss  per  round-trip  transit.  The  remaining  terms  represent  the 
traveling  wave  sources  due  to  the  aperture.  The  w  results  from  a  time- 


derivative  in  the  aperture  current  dipole  moment  representation,  with 
and  ag  the  required  components  of  the  dyadic  aperture  polarizabilit¬ 
ies.  The  I  h  >  and  I  e  >  vectors  are  associated  with  the  multiconduct- 

‘  x  y 

or  transmission  line  (MTL)  fields  at  the  aperture  as  required  for  the 
source  determination  using  the  reciprocity  theorem. 

For  a  particular  problem  as  shown  in  Figure  12  of  Chapter  IV  the 
radian  frequency  w  is  assumed  known  and  the  loss  must  be  estimated. 
This  loss  is  due  to  power  absorbed  in  the  terminations,  or  line,  and 
radiated  in  other  neglected  modes  or  back  through  the  aperture.  For 
simplicity,  we  shall  set  equal  to  only  the  line  loss  which  will  be 

estimated  at  ldB/lOOft. ,  typical  of  standard  transmission  lines  in  the 
high  frequency  (HF)  range. 

In  order  to  facilitate  the  field  calculations  associated  with  the 

sources,  the  medium  is  assumed  to  be  homogeneous.  In  general,  we  shall 

use  the  parameters  of  a  vacuum,  and  which  provide  a  maximum 

0 

value  of  v^  to  be  3  *  10  m/s.  The  bounds  on  the  polarizabilities  are 
obtained  from  the  geometry  of  the  aperture.  In  general,  these  bounds 
may  be  determined  from  the  polarizabilities  of  an  ellipse  which  cir- 

SC”  SC” 

cumscribes  the  aperture.  The  and  are  obtained  from  the 

exterior  problem  which  is  assumed  to  have  been  solved. 

The  remaining  quantities  A . ,  |  h  >,  and  |e  >  are  determined  from 

lx  y 

the  geometry  of  the  MTL.  The  equations  defining  the  MTL  are  given  in 


Appendix  D  as 


c  is 


The  parameters  are  the  eigenvalues  of  L  .  If  the  medium  is  a 
vacuum  as  assumed,  these  are  simply  the  eigenvalues  of  c^K  where 

g 

the  speed  of  light  3  *  10  m/s.  The  induction  coefficient  matrix  K 
may  be  determined  by  either  solving  n  boundary  value  problems  using 
Laplace's  equation  for  the  n  wires  or  determining  the  potentials  due 
to  charge  distributions  on  the  wires.  The  latter  has  been  used  where 
a  line  charge  has  been  assumed  and  the  associated  potentials  been 
determined.  These  results  are  used  to  determine  the  matrix  K  from  the 
vector  equation 


V  >  =  P  |  Q  > 


(12) 


with  the  elements  of  |  V  >  and  ]  Q  >  associated  with  the  corresponding 
wires.  The  matrix  K  is  the  inverse  of  P  from  which  the  ^  may  be  found. 

In  determining  _P,  the  electric  field  in  the  medium  is  required 
which  is  also  used  to  explicitly  determine  at  the  aperture  for  each 
Q^.  The  corresponding  magnetic  field  component  Hx  may  simply  be 
obtained  as  (-E^/nQ),  where  nQ  is  the  characteristic  impedance  of  a 
vacuum  given  by  120ir. 

To  complete  the  problem,  modes  on  the  MTL  must  be  defined.  In 

2 

general,  the  modes  are  associated  with  the  eigenfunctions  of  (-u>  LK). 

-1  2 

However,  in  a  homogeneous  medium  L  =  c  K  and  we  need  the  eigen- 
2 

vectors  of  ( - ^  I)  where  I_  is  the  identity  matrix.  In  this  case,  we 

c 

are  free  to  choose  the  eigenvectors.  It  is  standard  to  choose  one  of 


the  voltage  eigenvectors  for  an  n-wire  line  as 
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bl> 


1 


(122) 


which  corresponds  to  the  bulk  mode,  and  the  other  nodes  |  >  as  ortho¬ 

normal  differential  nodes.  The  corresponding  charge  distributions 
|  >  are  given  by 

|  Qi  >  =  K  ]  b±  >  (123) 

which  results  in  aperture  field  E  ^  for  each  |  >  .  The  components  of 

|  e^.  >  and  |  h^  >  are  normalized  for  unit  power  flow.  Thus 


/  2  E  . 
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(124) 


(125) 


and 

h  >  - - |  e  >  ,  (126) 

x  n0  y 

With  these  vectors  and  equation  (66) ,  we  are  prepared  to  determine  an 
upper  bound  on  the  frequency  domain  voltage  and  current. 

Alternately,  the  |  b^ >  may  be  chosen  as  the  unit  normalized  eigen¬ 
vectors  of  L  \  then  (123)  becomes 


c 


(127) 
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The  denominator  of  (124)  similarly  becomes  /xi/v  .  Since  E  .  is  pro- 

xi  71 

portional  to  — r,  (124)  can  be  simplified  to 
c^ 


where  the  prime  denotes  the  electric  field  due  to  a  charge  vector  |  b^  >. 

The  first  problem  to  consider  in  this  chapter  consists  of  two 
parallel  thin  wires  in  a  vacuum  as  shown  in  Figure  23.  A  particular 
frequency  is  not  chosen,  but  rather  the  answer  is  given  as  proportional 
to  a).  Using  the  suggested  ldB/lOOft  estimate  for  the  loss  along  with 
lengths  l ^  and  l ^  of  5  and  7  meters  respectively,  (1-e  °T)  ^  becomes 
(11.539).  These  are  the  same  lengths  as  the  example  of  Chapter  VI. 

Also  assuming  the  same  circular  aperture,  and  ag  are  given  by 
(1.333  x  10  6)  anci  (6.666  x  10  ^ )  respectively. 

To  determine  the  eigenvalues  and  field  vectors,  we  consider  the 
approximate  line  charge  equivalent  configuration  of  Figure  24.  To 
compute  the  P  of  Equation  (121) ,  we  may  sequentially  allow  only 
charge  and  charge  to  be  nonzero.  From  each  of  these  line  charges, 
the  electric  fields  may  be  determined.  We  obtain  the  corresponding 
voltages  by  integrating  the  electric  field  from  the  ground  plane  to  the 
wire  edges  (dotted  line  about  line  charge)  along  both  paths  and  L^. 
From  the  latter,  the  matrix  elements  of  £  are  given  for  thin  wires  as 


Aperture 


Figure  23.  Cross-section  of  two  wires  behind  an  aperture. 


Figure  24.  Image  charge  equivalent  of  two  wires  above 


a  ground  plane. 
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and 
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Thus 
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(131) 


g 

where  c  is  the  speed  of  light  (3  *  10  m/s) . 

It  is  easily  shown  that  the  eigenvalues  of  L  ^  are  given  by 


*1  ’  c  /<PU+P12> 


(132) 


with  the  associated  eigenvectors  given  by 


b.  >  =  -±- 

1  n 
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+i 


(133) 


If  the  charges  on  the  two  wires  are  designated  as  and  C^,  then 
it  is  easily  shown  that  the  electric  field  at  an  aperture  symmetrically 
located  between  the  wires  is  given  by 
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Choosing  |  >  to  have  the  plus  sign  we  have 
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(135b) 


We  note  that  is  the  minimum  eigenvalue  required  in  Equation  (66) . 
Substituting  this  information  into  (66)  we  obtain 
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tt  (d  +a  ) 


+  !Eysc'l3. 


(136) 


If  a  =  d  =  1cm,  this  becomes 


|A  V|  <  9.794  xlO-4  k  [2 |q  H  Sc"|  +  |E,SC“|l 
1  'max  ol  1  o  x  '  1  y  1 


where  k  is  the  wave  number  w/y  e  .  Using  the  H  SC  and  E  ai"  of  (118) 
o  oo  x  y 


sc- 


with  the  associated  data,  this  bound  is 
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AV |  <  334  k  V 

'max  o 

or  equivalently 

I AV|  <  oo(l.ll  xio"6)  V  .  (137) 

1  max 

This  is  approximately  five  times  the  level  of  the  parallel  plate 
problem  of  Chapter  VI  with  a  2  cm  spacing. 

To  obtain  the  time-domain  upper  bound  for  the  vacuum  filled  line 
we  note  that  the  factor  u>/(l-e  °^)  of  equation  (66)  must  simply  be 
replaced  by 

m 

I 

j=l 

to  obtain  equation  (108).  In  (108),  the  loss  has  been  neglected  in 
comparison  to  the  decay  of  the  short  circuit  fields  as  specified  by 
for  each  mode.  Using  the  data  of  Kajfez  in  equation  (120),  the  time 
domain  bound  becomes 

| AV |  <  10.98V.  (138) 

1  1  max 

The  second  problem  of  interest  is  shown  in  Figure  25  with  wire  2 

close  to  wire  1  and  much  smaller  with  both  radii  much  less  than  d.  This 

problem  models  a  wire  at  the  surface  of  a  wire  bundle  with  an  aperture 

centered  below.  If  one  assumes  that  the  only  change  from  the  previous 

problem  is  in  the  cross-sectional  geometry,  only  the  new  X^, 

je  >,  and  |  h  >  are  needed  for  equation  (66)  to  be  used.  Fundamental 
y  x 

to  these  quantities  is  the  determination  of  £  in  (121). 

To  determine  P  ,  a  charge  set  of  =  0  and  =  1  is  first  con¬ 
sidered.  Due  to  the  proximity  of  the  wires,  the  model  of  Figure  26  is 


Aperture 


Figure  25.  Cross-section  of  thick  and  thin  wires  behind 
an  aperture. 


Figure  26.  Image  charge  equivalent  for  only  small  wire 
charged. 
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required  with  the  center  charges  representing  the  cylindrical  images  of 
(?2  in  wire  1  and  the  left  charges  balancing  the  center  charges  to  give 
=  0.  Computing  the  resulting  voltages  and  between  the  wires 
and  the  ground  plane,  the  corresponding  forms  of  P ^  and  P22  are 
obtained  as 


12  4ire  a 
o 


“ 

2d  -rL 

1  , 

In 

+  4  In 

A 

_ 

2  ,  .2 
r^  +  r 


(2d  -r1)2+r’2 


4ire 


In 


(2d-ri)2+4a2 

2.,  2 

r^  +4a 


and 


22  8ire  a 
o 


In 


(2d  -r2)2  +4a2 


2  .  2 
r2  +  4a 


+  In 


r22  +  (2a-r')2 
(2d-r2)2+(2a-r’)2 


2ire 


In 


2d  -  r„ 


where  r^  and  r2  are  the  radii  of  wires  1  and  2  and  r  =  r^  /2a  is  the 
location  of  the  image  charge. 

For  P22  and  ^22’  no  *ma8e  is  needed  in  wire  1  and  the  line  charges 
are  simply  =  1  and  Q2  =  0  to  give 

2d  -  r, 


'  11  2ire 


and 


P„  = 


21  4tt£ 


In 


In 


(2d  -r2)2  +Aa2 

2  ,/  2 
r„  +  4a 


Dimensions  are  chosen  similar  to  the  previous  problem  as  d  =  1cm 


r^  =  0.5cm,  ~  0.1cm,  and  a  =  0.35cm.  The  image  position  r  is  given 

by  0.3571cm.  The  resulting  matrices  £  and  L  ^  are  given  by 


—  ~  2ire 


1.0986 

1.0521 


0.7821 

3.3940 


and 


L_1=  — 
~  yo 


1.1680 

-0.3621 


-0.2692 

0.3781 


The  required  eigenvalues  and  eigenvectors  of  L  ^  are  given  by 


.  _  1.6939  ,  _  8.0205 

A1  VJ  *  A2  “  y 


and 


0.2870 

-0.9275 

= 

,  |b2> 

0.9579 

L 

0.3738 

where  ]b^>  represents  the  bulk  mode. 

The  aperture  is  located  half  way  between  the  centers  of  the  wires 

as  shown  in  Figure  25.  The  aperture  electric  field  for  each  mode  is 

needed  for  |e  >  and  ]h  >.  For  mode  1,  the  charge  vector  is 
y  x 


iQ>  “  I_1  1^  >  =  2irec 


0.0773 


0.2583 


To  compute  at  the  aperture,  must  be  modeled  by  three  charges  above 
the  ground  plane  as  was  done  in  the  evaluation  of  and  ? 
in  mind  the  general  formula  for  is 


With  this 
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E  =  — 

y  ire 


r1^2  d 

(Ql+Q2+  -j^-)  2  2 

a  +d 


"  Qo  oT 


2  2a  r ,2  .  .  2,.  .2, 

[d  +  (a  -  rx  /2a)  ] 


to  give 


Ey  =  -55.78  V. 


Similarly  for  mode  2, 


|q  >  *  2we 
1  o 


-1.1839 

0.4772 


and 

Ey2  =  133.35  V. 

2 

From  equation  (128),  |e  >  becomes  with  the  (A./c  )  scale 

J 

-1176. 8~ 

1292.9 

Substituting  into  (66),  the  upper  bound  becomes 

|AV|  <  0.8025  toe  [2  i  n  H  SC_  |  +  |E  SC"|]. 

'  'max  o  1  o  x  '  1  y  1 

For  the  incident  form  of  Kajfez,  this  bound  is 

I AV I  <  u(2.41  x 10-6)  V. 

1  'max 

In  the  time  domain  this  bound  becomes 

|AV|  23.84  V. 

1  'max 

The  substantial  increase  in  these  bounds  compared  to  the  parallel  plate 
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problem  is  seen  as  a  primary  result  of  the  five-fold  increase  in  the 
L2  norm  of  |e^>.  This  is  a  consequence  of  the  closeness  of  the  wires 
and  aperture  and  the  stronger  interaction. 

The  last  problem  is  depicted  in  Figure  27.  The  wire  has  been 
placed  at  the  center  of  the  structure  and  has  a  radius  equal  to  b/100. 
The  wire  is  labeled  as  conductor  1  with  the  upper  plane  as  conductor 
2.  If  =  0  and  Q2  =  1,  the  wire  may  be  neglected  in  computing  P  2 
and  P2?  to  obtain 


12  e  L 
o 


22  e  L  * 
o 


The  computation  of  and  P21  is  more  complicated  since  the  images  of 
the  wire  in  the  surfaces  must  be  included.  The  values  P^  and  P21  may 
be  obtained  by  subtracting  times  P^2  °r  P22  ^rom  t*ie  potentials  for 
a  wire  between  two  grounded  planes.  This  removes  the  effect  of  charge 
on  conductor  2.  For  the  grounded  planes  case  with  the  wire  charge 
equal  to  1  C/m,  the  potential  between  planes  is  zero  and  the  wire 
potential  is  obtained  from  an  infinite  sum  including  the  images.  If  only 


the  lower  image  is  used,  the  wire  potential  is  approximated  by 
(4. 595/2TreQ) .  The  additional  images  have  a  minor  contribution  changing 
this  value  to  (4. 0466/2ireo) .  Combining  these  results  for  the  geometry 
given  with  L  *  10b,  one  obtains 

P 0.6854  0.05" 

P  =  -~- 

Zo  0.05  0.1 


life 


The  corresponding  eigenvalues  and  eigenvectors  are 


,  1.45  ,  10.44 

X1  IT  *  X2  

o  o 

and 


0.9963 

-.0845 

= 

.  |b2>  = 

0.0845 

0.9964 

The  aperture  is  located  directly  below  the  wire  for  maximum 


coupling  and  d  is  set  equal  to  1cm.  The  image  summation  for  E^  can  be 


computed  exactly  for  this  problem  to  give  the  electric  field  at  the 
aperture  as 

i—  —i 

Qo 


E  =  Q1 
y  X1 


2Le 


4de 


Lc 


Using  the  |b^>  as  charge  distributions,  the  primed  fields  of  (128)  are 
obtained  as 


E  '  =  -22.84/e 
y±  o 


and 


h2  =  "3-°8/£o 

which  give  an  |e  >  of 


Since  |e  >,  jh  >,  and  A  are  the  only  parts  of  (66)  to  be 
y  x  1 

modified  from  the  last  problem,  the  upper  bound  is  obtained  as 

1 AV |  <  u(1.197  x 10-6)  V 

1  1  max 

with  an  equivalent  time  domain  bound  of 

AV  j  <  11.84  V. 

'max 

The  primary  observation  of  interest  which  seems  to  link  these  three 
problems  is  that  for  similar  dimensional  relationships,  the  bound  takes 
on  the  same  order  of  magnitude.  This  would  suggest  the  applicability 
of  canonical  problems  as  bounds  for  more  complex  problems. 
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CHAPTER  VIII 
CONCLUSIONS 

This  report  developed  a  computational  technique  for  upper  bounds 
on  signal  levels  at  terminations  of  multiconductor  transmission  lines 
(MTL)  located  behind  an  aperture-perforated  conducting  screen.  This 
was  accomplished  in  several  stages.  First,  the  electromagnetic  coupling 
through  small  aperture  was  described  using  the  concept  of  aperture 
polarizabilities.  The  idea  of  coupling  was  then  extended  to  the  MTL 
excitation  behind  an  aperture  where  the  aperture  was  represented  by  the 
equivalent  current  dipole  moments.  A  source  of  traveling  waves  was 
Introduced  to  replace  the  aperture  and  the  aperture  coupled  energy. 

These  traveling  waves  transferred  the  energy  from  the  aperture  region 
to  the  terminations.  The  amount  of  transferred  energy  in  terms  of  the 
traveling  waves  was  found  by  introducing  a  signal  flow  graph  of  the 
whole  system.  It  was  assumed  that  the  aperture  had  no  other  signifi¬ 
cant  interactions  with  the  MTL  as  might  be  described  by  additional 
reflections  in  the  aperture  region. 

Having  computed  the  amount  of  energy  at  the  terminations,  a  trans¬ 
formation  was  used  to  represent  the  traveling  waves  in  terms  of  the 
voltages  and  currents.  The  mathematical  properties  of  vector  norms  and 
their  associated  matrix  norms  were  discussed.  The  idea  of  two-norm  was 
used  to  formulate  upper  bounds  on  the  voltages  ar.d  currents  at  the 
terminations.  At  this  stage,  it  was  assumed  that  the  terminals  were 
passive  and  the  losses  associated  with  each  mode  traveling  on  the  MTL 
were  of  exponential  form.  Although  the  radiation  and  termination  losses 
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were  not  explicitly  taken  into  account,  it  was  mentioned  that  these 
losses,  if  non-negligible,  could  be  considered  by  modifying  the  term 
due  to  propagation  loss.  Further,  after  some  mathematical  manipulations, 
it  was  found  that  the  bounds  depend  primarily  on  the  source  and  termi¬ 
nation  local  geometries.  Finally,  considering  the  fact  that  the 
maximum  voltage  between  any  two  lines  might  be  the  sum  of  the  indivi¬ 
dual  voltages  to  ground,  an  additional  factor  of  two  was  Incorporated  in 
the  upper  bound.  The  upper  bound  on  the  termination  voltages  was  formu¬ 
lated  in  the  frequency  domain.  An  upperbound  on  the  currents  in 
frequency  domain  was  simply  obtained  from  this  voltage  bound. 

This  idea  was  extended  to  the  time  domain  where  it  was  assumed  that 
the  waveforms  launched  from  the  sources  were  non-dispersive  (quasi-TEM) , 
thus  reaching  the  terminations  without  any  distortion.  Several 
assumptions  were  made  regarding  the  medium,  the  termination,  and  the 
modal  propagation.  It  was  assumed  that  the  elements  of  the  termination 
are  all  real  and  constant  and  the  propagation  constant  of  all  the  modes 
are  equal  to  that  in  a  vacuum.  The  multiple  reflection  phenomenon  was 
taken  into  account  and  its  formulation  was  simplified  by  using  the  pole 
expansion  of  the  incident  field  characterized  in  the  Singularity 
Expansion  Method  (SEM) .  Using  vector  and  matrix  norms,  an  upper  bound 
for  the  voltage  and  the  current  was  formulated.  It  was  observed  that 
Che  upperbounds  basically  depend  on  the  geometry  of  the  system,  the 
properties  of  the  external  time-varying  fields,  and  the  shape  of  the 
aperture. 

Having  formulated  the  upperbounds  on  the  signal  levels  in  both  time 


and  frequency  domain,  a  comprehensive  example  was  presented  in  order 
to  elaborate  the  use  of  the  technique.  The  parameters  and  the  geometry 
were  chosen  exactly  the  same  as  a  problem  introduced  by  Kajfez  [12]  in 
order  to  facilitate  the  comparison  of  results.  It  was  found  that  the 
maximum  voltage  in  Kaj fez's  example  for  a  moderate  mismatched  termina¬ 
tion  was  found  by  a  multiplicative  factor  of  approximately  ten.  It  was 
noted  that  for  upperbound  calculations  under  the  worst  possible  con¬ 
ditions  a  factor  of  ten  was  in  reasonable  agreement  with  his  results. 

The  reason  being  that  at  several  stages  of  computation  various  types  of 
inequalities  were  used  that  doubled  the  exact  value.  Also,  the  velocity 
of  all  the  modes  were  assumed  to  be  the  velocity  of  light  in  a  vacuum 
which  were  greater  than  the  actual  velocity  of  the  modes  in  their  res¬ 
pective  media,  contributing  to  the  increase  of  the  voltage  at  the 
termination.  In  addition,  a  factor  of  two  was  used  to  account  for  the 
potential  difference  between  any  two  wires  and  not  just  to  ground.  In 
the  whole,  It  was  observed  that  final  results  obtained  by  using  the 
developed  techniques  were  in  good  and  reasonable  agreement  with  Kaj fez's 
results. 

Several  other  problems  were  also  approached  directly  by  the  techniques 
developed.  Two  bundles  of  wires  over  a  ground  plane  were  bound  above 
by  modeling  them  as  two  wires  over  a  ground  plane.  A  wire  at  the  surface 
of  a  bundle  of  wires  was  modeled  by  separate  thick  and  thin  wires  over  a 
ground  plane.  A  last  example  considered  a  wire  between  parallel  plates. 

The  problems  all  had  comparable  dimensions  but  interestingly  had  similar 


voltage  bounds. 
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There  are  several  suggestions  for  future  research.  It  would  be 
desirable  to  find  signal  levels  on  wires  behind  composite  panels 
which  couple  the  external  field  by  a  diffusion  mechanism.  Also  of 
some  importance  is  the  problem  of  wires  passing  through  an  aperture 
corresponding  to  antennas  or  aircraft  control  cables.  Another  area 
of  research  would  be  large  aperture  problems  such  as  conformed 
antennas  or  large  windows. 


Appendix  (A) 

POLARIZABILITY  OF  SMALL  APERTURES 

The  diffracted  field  in  the  vicinity  of  an  aperture  depends  on  the 
excitation  field  and  upon  the  shape  and  size  of  the  aperture.  The 
moments  of  the  equivalent  dipoles  are  related  to  the  components  of  the 
known  exciting  field  through  special  constants  of  proportionality 
called  the  aperture  electric  and  magnetic  polarizabilities. 

Due  to  the  importance  of  the  use  of  dipole  moments,  a  great  deal 
of  attention  has  been  devoted  to  determination  of  the  polarizabilities. 
Cohn  in  [AlJ  and  [A2]  has  experimentally  determined  the  constants  for 
several  shapes,  while  Van  Bladel  [A3],  has  computed  by  numerical 
methods  the  polarizabilities  for  a  rectangle,  diamond,  cross,  and  a 
rounded-off  rectangle.  Table  (Al)  gives  the  electric  and  magnetic 
polarizabilities  for  a  circle  of  radius  (R) ,  an  ellipse  of  eccentricity 
(e) ,  and  a  narrow  ellipse.  The  following  remarks  should  be  made 
concerning  the  cable: 


1) 


a 

m 


a 

m,  xx 


xx 


a 

m,yy 


3)  K  and  E  are  the  complete  elliptic  integrals  of  the  first  and 
second  kind,  as  in  [A5J. 

By  the  3tudy  of  the  data  given  by  De  Meulenaere  [A3],  when  an  in¬ 
accuracy  in  the  polarizabilities  of  10%  can  be  tolerated,  one  can 

3/ -» 

use  a  normalization  factor  of  [(aperture  area)  "],  and  use  elliptic 


polarizabilities  in  order  to  calculate  the  polarizabilities  of  the 
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rectangle  and  rounded-off  rectangle. 

The  polarizabilities  of  other  shapes  nay  be  obtained  by  either 
measurement  or  numerical  solution  of  the  quasi-static  aperture  problem. 


Table  (Al):  Aperture  Polarizabilities  [13) 


major  diameter  of  ellipse 


Appendix  (B) 

APERTURE  REPRESENTATION  BY  A  PAIR 
OF  CURRENT  DIPOLES 

Consider  Figure  (Bl)  where  an  aperture-perforated  screen  is 
shown  [Bl].  In  order  to  evaluate  the  scattered  field  in  the  interior 
region,  one  only  needs  to  know  about  the  tangential  electric  field 
£c  over  the  aperture-  It  is  convenient  to  close  the  aperture  with  a 
metalic  lid  and  place  a  magnetic  surface  current  density  J ^  over  it 
given  by 

Jsm  -  Ec  x  n  (Bl) 

where  n  is  the  normal  vector  as  shown  in  Figure  (B2) .  The  vector  E  is 
the  total  aperture  field  typically  obtained  by  solution  of  a  quasi¬ 
static  integral  equation  [B2]. 

Now,  one  can  invoke  the  reciprocity  theorem  to  compute  the 
scattered  field  at  any  point  in  the  interior  region  due  to  surface 
magnetic  current.  For  the  purpose  of  this  discussion,  the  reciprocity 
theorem  as  introduced  in  [B3]  will  be  used.  As  shown  in  Figure  (33), 
the  scattered  field  f(r')  at  a  point  (r ' )  is  to  be  computed  with 
respect  to  the  origin  which  is  taken  at  the  center  of  the  closed 
aperture  region.  By  introducing  a  unit  magnitude  electric  dipole  P 
at  the  point  (r' )  [Bl],  as  3hown  in  Figure  (B3) ,  one  can  write  the 
reciprocity  theorem  as 

f(r').  a  -  (J  .  E,  -  J  .  il  )  dV  (32) 

p  e  o  m  b 

V 

where  ?  is  an  electric  current  dipole  given  by 

P  -  a  5  (r  -  r’)  (33) 

o 
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Figure  (Bl) .  Aperture-perforated  screen  [8]. 


%  -  e 

E_xn  *  J  „  _  «<■  —  -  . 

c  3  ©  ®  ® _ 

|  Oucside  Region 

Figure  (B2) .  Aperture  replaced  by  j”  [8], 
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Figure  (B3) .  Evaluation  of  scattered  field  [8]. 


and  ap  is  a  unit  vector  characterizing  the  orientation  of  the  testing 
dipole  r.  The  vectors  and  are  the  fields  produced  fay  the  source 
P,  evaluated  at  the  surface  of  the  closed  aperture,  and  and  are 
volume  electric  and  magnetic  current  densities  which  are  the  sources 
of  the  field  scattered  by  the  aperture.  In  the  case  of  this  discussion, 

«*•  -►  Qj 

J  is  equal  to  zero  and  equation  (B2)  for  the  surface  current  J  becomes 

► 

E(r')  .  tp  -  -  (Js“  .  H^  dS.  (B4) 

S 

From  Maxwell's  equations,  one  can  write 

Vx^  -  jiue  ?  (B5) 

The  vector  can  be  expanded  in  Taylor's  series  as 

\(t)  -  1^(3)  +  v  .  (VH)  -  +  ...  (B6) 

where  the  higher  order  terms  are  assumed  to  be  negligible.  One  can 
write  r  .  (7H) | -  as 

r  •  (VH) | -  -  (7H) | -  .  r  -  rx  (7xH)|-.  (B7) 

The  term  (7H)|-  •  r  gives  rise  to  a  quadrupole which  is  neglected  [32]. 
Substituting  equation  (B5)  in  equation  (B7) ,  one  may  rewrite  equation 
(B6)  as 

H^Cr)  -  H^(o)  -  r  x  [jwe  £^(3)]  .  (B8) 

If  equation  (B8)  i3  substituted  in  equation  (B4) ,  one  obtains 

E(r')  •  ap  -  -[H^(o)-  jucr  X  Eb (3) ]  •  J,md S 

Js  (B9) 

One  mav  define  C  and  C  as  electric  and  magnetic  diDole  current 
e  m  ‘ 


moments  respectively  by 
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C  *  iue 
e 

S 

C  -  j“  ds  .  (Bll) 

m  o 

s 

With  these  definitions,  equation  (B9)  can  be  written  as 

E(r')  .  ap  -  -  ^(o)  .  Ce  -  ^(o)  .  Cm  (B12) 
By  comparison  of  equations  (B2)  and  (B12) ,  one  may  write  equivalent 
currents  as 

-  <5(y)  JmS  -  C  d(r)  (B13) 

m  m  m 

Je  -  5 (y)  JeS  -  C06(r)  (B14) 

where  <5(r)  is  the  three  dimensional  Dirac  delta  distribution.  The 
equivalent  dipole  current  moment  representation  of  an  aperture  is 
shown  in  Figure  (B4) . 


Figure  (34).  replaced  by  a  pair  of  current  dipoles. 


Appendix  (C) 

DIRAC'S  NOTATION  OF  MATRICES 

The  following  material  is  an  introduction  to  matrix  algebra  using 
Dirac's  Notation  which  was  introduced  in  Quantum  theory  by  P.A.M.  Dirac. 
The  presented  material  is  a  summary  of  B.  Friedman  [Cl],  S.  Gasiorowicz 
[C2],  and  A.  Messiah  [C3]. 

In  Dirac's  Notation,  an  N-dimensional  column  vector  is  denoted  by 


When  matrix  A,  an  N x  N  matrix  operates  on  vector  |x>  ,  the  resulting 
operation  is  another  column  vector  [V>.  By  elementary  properties  of 
matrix  multiplication,  the  kth  element  of  the  vector  |V>is  computed  as 

<C1) 

The  complex  conjugate  transpose  of  a  column  vector  |x>  is  a  row 
vector  denoted  by  <  x |  whose  elements  are  the  complex  conjugate 
elements  of  jx>  ordered  as 

<  x  |  -  (x*  x2  ...  Xjj  )  (C2) 


where  (*)  means  complex  conjugate. 


The  (t)  nocation  is  used  for  complex  conjugate  transpose  of  A 


given  by 

A+  - 


Using  matrix  properties  for  conjugation  and  multiplication  the 
following  product  is  . obtained: 

<  M  |  *(AB  |x>  y  »  <  x  |  Bf  A1-  (C3) 

where  <M  |  is  a  Nxl  row  vector  and  B*  and  A1’  are  both  NxN  matrices. 

One  of  the  most  important  types  of  matrices  is  the  Hermitian 
matrix  which  is  equal  to  its  own  complex  conjugate  transpose,  i.e., 

Ar  -  A 

A  Hermitian  matrix  has  the  following  properties: 

2 

1)  Instead  of  (N  )  distinct  elements,  it  only  has  *iN(N+l)  distinct 
elements. 

2)  The  elements  on  the  main  diagonal  are  real. 

3)  The  elements  which  are  located  symmetrically  across  the  main 
diagonal  are  complex  conjugate  of  each  other. 

Multiplication  of  a  row  vector  with  a  compatible  column  vector 
results  in  a  complex  number  as 


*  * 


'll 

*111 

* 

* 

12 

*N2 

'* 

111 
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10 


<  x  !  y> 


?N 


(C4) 


It  may  be  observed  that  <  x|y  >  «<y|x  >*  .  If  < V  ]  ■  <  x|A‘f  then 
<  V|y  >  ■  <  xjA^|y> ,  and  the  result  of  such  an  operation  is  scalar. 

Sometimes  it  may  happen  that  the  order  of  multiplication  in  (C4)  is 
interchanged,  then  the  result  is  a  N*N  square  matrix 


ly><  x| 


(C5) 


One  of  the  benefits  of  representing  matrices  in  Dirac's  Notation  is 
that  by  observing  the  position  of  symbols,  one  can  readily  determine  the 
format  of  the  resulting  operation.  For  example,  assuming  that  3  is  a 
scalar,  |x  >  and  |y  >  are  vectors,  and  A  is  a  square  matrix,  then 
<  x|A|y  >  and  S  <  y|x  >  are  scalars,  a[x  >  is  a  vector,  and  jx  ><  y| 
is  a  square  matrix. 


Appendix  (D) 


MULTICONDUCTOR-LINE  FORMULATION 

The  computation  of  voltages  and  currents  on  a  multiconductor 
transmission  line  (MTL)  in  terms  of  eigenvectors  have  been  analyzed  by 
Amemiya  [Dl]  and  Marx  [D2J ,  This  approach  has  been  used  for  both 
transient  and  steady-state  waveforms.  D.  Kajfez  [D3]  has  extended  the 
eigenvector  treatment  by  using  a  simultaneous  diagonalization  of  two 
matrices.  This  section  presents  the  main  ideas  of  method  of  diagonal¬ 
ization  by  using  only  one  composite  matrix  instead  of  two  separate 
matrices. 

Consider  the  N-conductor  transmission  line  with  a  conducting 
ground  plane  as  shown  in  Figure  (Dl) .  The  reference  direction  has 
been  chosen  such  that  real,  positive  values  of  and  1^  represent 
power  flow  in  the  positive  z  direction. 

References  [Dl],  [D2J,  and  [D4]  derive  the  following  two  formulas 
for  a  lossless  MTL  in  sinusoidal  steady-state  analysis  with  isotropic, 
nonmagnetic  dielectrics: 


_d_ 

dz 

|V  >  -  -jwL|  I> 

(Dl) 

_d_ 

dz 

|l  >  -  -ju)K|  V> 

(D2) 

The  matrix  is  called  the  induction  matrix  and  the  matrix  K  is  the 
induction  coefficient  matrix. 

There  are  two  types  of  propagation  on  these  lines.  If  the 
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dleleccric  is  homogeneous,  transverse  electromagnetic  (.TEM)  modes 
propagate.  If  the  dielectric  is  inhomogeneous,  hybrid  electromagnetic 
(HEM)  modes  propagate  [D5].  In  HEM  mode  both  the  electric  and  magnetic 
vectors  have  a  longitudinal  component.  At  low  frequencies,  the  longi¬ 
tudinal  component  is  small  and  may  be  neglected,  giving  rise  to  quasi- 
TEM  modes. 

Using  the  analysis  presented  by  Friedman  [D6],  the  first  step  is- 
to  decouple  equations  (Dl)  and  (D2)  by  taking  derivatives  of  both  sides 
and  substituting  in  the  same  equations  to  obtain 

j2  , 

-=-=■  |V  >  -  -tTLK  |  V>  (D3) 

dzZ 

A  , 

-2-r  |l  >  -  VKL  |  I>  .  (D4) 

d  tT 

As  and  K  represent  stored  energy  of  a  passive  network,  they  are  both 

positive  definite  and  thus  qualify  for  the  method  of  Friedman  [D6]. 

This  method  of  solution  for  |v>  and  j  I>  consists  first  of  obtaining  a 

set  of  orthonormal  eigenvectors  |  ?i>  and  corresponding  eigenvalues 
2  2 

(  -  )  for  the  matrix  LjC  .  Since  JL  and  K  are  positive  definite, 

2 

LK  is  also  positive  definite  and  the  6^  '  s  are  all  positive.  The 
eigenvalue  problem  is  given  by 

-<Ak  |ci>-  -Si2ki>  •  (D5) 

The  eigenvalues  may  be  obtained  from 

det[-w2L  K.  +Si2  IJ  -  0  (D6) 


where 


•i} 


CD7) 


and  <5^  is  Che  Kroneker  delta. 

One  may  now  form  a  square  matrix  G  whose  columns  are  the  | > 

£  -  r  |?1>,  |52> .  UN>]  (D8) 

This  matrix  may  be  used  to  diagonalize  _B  as 

A  5  G+ B  G  -  Diag[  -  6* . -  SN2  ]  (D9) 

where 

GGf  -  I  .  (DIO) 

The  matrix  G  is  an  orthogonal  real  matrix,  (.+)  denotes  the  transpose  of 
the  matrix,  and  _I  is  the  identity  matrix. 

The  corresponding  eigenvalue  problem  of  the  decomposed  matrix  A 
becomes 

,2 

|Y  >  -  A  |  Y>  (Dll) 

dz 

where 

|Y  >  -  G+|V  >  (D12) 


3ecause  of  the  fact  that  A  is  real  and  diagonal,  equation  (Dll)  can  be 
decomposed  into  a  second  order  differential  equation  for  each  mode  as 


?i 


(D13) 


The  solution  is  of  the  exponential  form 
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■  y.,+  e  'i  +  y^~  e 


+jSiz 


(D14) 


where  S.  *  — - ■  ,  and  v.  is  the  velocity  of  the  ith  mode.  This  velocity 
i 

relation  can  be  shown  by  the  usual  method  of  multiplying  (D14)  with 
e+^ut  and  then  taking  the  real  part  in  order  to  obtain  the  rate  of 
change  of  the  phase.  In  equation  QD14),  and  y^  refer  to  the 
amplitudes  of  positive  and  negative  travelling  waves  in  the  z  direction. 
|Y  >  can  now  be  written  as  (D15) 


|Y>  -  f  y  |u  > 
i-1 


(D15) 


where  |u^>  denotes  a  column  vector  having  all  zero  elements  except  for 
the  ith  element  which  is  unity. 

Equation  (D12)  can  now  be  used  to  find  jV>  as 


IV  > 


N 

I  Yii?i> 


(D16) 


i-1 


Similarly,  |l>  can  be  obtained  from  equations  (Dl)  and  (D16)  as 


2  -i  dvj  -i  , 

ir>- i‘i> 


(017) 


The  above  equation  can  also  be  written  as 

|I  >  -  !  C  V  1  I='l+  e'J8lZ  -  r[  e+38l']l'1Ui> 

i-1  1  X 


(D18) 


From  these  expressions  for  voltage  and  current,  the  total  power 
carried  on  a  MTL  is  the  average  real  power  defined  by 
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P  -  h  Re  <  V 1 1  > 


©19) 


Substituting  equations  ©18)  and  (D16)  into  (D19),  P  becomes 


P  -  h  Re 


li,j-i  “  1  j  1  j  >  • 


(D20) 


where  C*)  denotes  the  complex  conjugate.  If  y^"  is  equal  to  zero  then 
the  power  travelling  in  positive  z  direction  is 

S  s, 


(D21) 


It  is  useful  to  relate  the  positive  traveling  power  (P  )  to  a 
power  wave  (a^)  such  that 


,  « 
?+  •hi 


(D22) 


J-l 


assuming  that  the  travelling  modes  are  orthogonal.  In  this  case,  ai  Is 


given  by 


li  *  , -  \/<  I  ^i  I  **i  ‘ 


(D23) 


where  the  X^'s  are  the  eigenvalues  of  the  matrix  L.  Generalizing,  the 
total  power  can  be  written  as 
N 

I 

i-1 


P  *  4  I  C  i  ai|2+  |bi|2  ] 

j _ s  ^  ^ 


(D24) 


where  the  b^'s  are  defined  by 


The  voltage  vector  |V>  may  be  written  in  terms  of  the  voltage 


eigenvectors  j<^>,  normalized  to  unity  power  as 
N 

|V  >  ■  l  (a  +  b  e'*eiZ)i$  > 


(D26 


where 


.  > 


/  v 
v  i 


i  j -  ”i 


> 


(D27 


Similarly,  is  the  corresponding  normalized  current  eigenvector 

giving  rise  to 

N 

l 

i-1 

where 


I>-  f  -  b1^Si‘)|»i> 


(D28) 


(D29) 


It  can  be  easily  shown  that 


<  ♦Jfj  >  - 


(D30) 


In  order  to  transform  from  the  power  wave  formulation  back  to  the 
desired  voltages  and  currents,  one  defines  two  matrices  whose  columns 
are  the  voltage  and  current  eigenvectors  given  by 

Hy  -  (  [*N>  )  (D31 


Ill 


Mj  -  C  |  *1 > , • • • ,  |  > ) 

From  equation  (D30),  it  is  easily  shown  that 

-I  ^7  "  - 

and 

^  hT  - 1  • 


(D32) 


(D33) 


(D34) 


In  the  bounding  process  of  Chapter  IV,  an  inhomogeneous  medium 
would  be  replaced  by  a  homogeneous  medium  with  a  dielectric  constant 
less  than  or  equal  to  the  minimum  dielectric  constant  of  the  inhomo¬ 
geneous  medium.  This,  typically  leads  to  the  use  of  a  vacuum  model 
2  2 

with  the  (g ?"  »  -^-y)  replaced  by  *  ^y) ,  where  c  is  the  speed  of 

1  Vj^ 

light  in  vacuum.  In  this  instance,  equation  (D13)  becomes 


yi 


(035) 


The  corresponding  and  would  also  change  to  account  for  the 


homogeneity. 


-•SmJ-. 
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APPENDIX  (E) 


TRAVELING  WAVE  FORMULATION 


The  voltages  and  currents  on  a  MTL  as  described  in  Appendix  (D) 


»  I  -J6.z 

|v(z)>  -  l  e  +  bi  e 


+jSiz 


|I(.)>-  I  ...  1  -b,.  1  V  (E2) 

i-1  - 

where  1 4>i>  and  ! <i'i>  are  normalized  voltage  and  current  eigenvectors, 
a.^  and  b^  represent  the  amplitudes  of  the  ith  mode  waves  traveling  in 
the  positive  and  negative  (z)  directions  respectively,  and  (6^)  is  the 
propogation  constant  of  the  ith  mode. 

Using  My  and  matrices  as  introduced  In  Appendix  (D) ,  equations 
(El)  and  (E2)  can  be  written  in  more  compact  form  as 


I V(z)>  -  My  (F  (z) | a>  +  F(z) | b> ] 


i  I (z) >  »  Mt  [F  ( z) | a>  -  F(z)jb>] 


where 


F(z)  »  diag 


f  j  8j_z 

e  ,  .  .  .  , 


Equations  (E3)  and  (E4)  can  be  written  as 


F  (z) ja>  -  f  jMI|v(z)>  +  Myil(z): 


I(z) i b>  •  4  I M^ j V(z) >  -  My i I(z) > | 
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In  the  above  formulas,  (a^)  is  the  complex  amplitude  of  ith 
mode  such  that  the  power  contained  in  the  positive  (z)  direction  is 


(E8) 


Hence,  the  entire  power  transmitted  along  the  positive  (z)  direction 
is 


PlUal  ‘  I  <a|a> 


(E9) 


Similarly,  the  total  power  contained  in  the  negative  (z)  direction  is 


P,  .  *  $  <b  b> 
Total  2  1 


(E10) 


where  it  is  assumed  chat  there  is  no  transfer  of  power  among  different 
modes . 

It  is  convenient  to  introduce  a  composite  MTL,  as  shown  in 
Figure  (El) ,  with  Z^  and  impedance  matrices  respectively  located  at 
z  =  Z,  and  z  *  ~l~  with  respect  to  the  origin  where  the  source  is 

4  J 

coupled  to  the  MTL  [El] .  The  amplitudes  of  the  ith  mode  introduced  by 

the  source  are  a  .  and  b  ..  The  total  amplitude  a.  of  the  ith  mode 
si  si  x 

consists  of  ag^  due  to  the  source  and  the  reflected  bg^  coming  from 
the  left  hand  termination  at  z  3  -Z^.  From  equations  (E6)  and  (E7) 
the  following  relations  can  be  obtained: 


!a<v>  *  I  |^iiv(^)>  + 

!b a4)>  -  \  |m^Iv(^)>  -  sjJ|i«4)>J 


(Ell) 


(E12) 
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By  Ohm's  law,  one  obtains 

|v(*4)>  "  >  I^4>>  <E13) 

Combining  equations  (£11) ,  (£12) ,  and  (£13) ,  one  obtains 

|b<*4)>  -  F4|a(£4)>  (£14) 

where  F4  is  called  the  reflection  matrix  given  by 

i4  "  i£h&i  *  i)  (^?4Mi  +  I)'1  (£15) 

The  matrix  F4  may  contain  non-zero  off-diagonal  elements  creating 
mode  coupling.  Similarly,  at  z  ■  -ty  one  can  obtain 

|a(-t3)>  -  £3|b(-£3)>  (E16) 

where 

£3  -  (mJZjMj  -  i) (M^Z3Mt  +  i)'1  (E17) 

One  can  also  look  at  the  problem  in  time-domain  assuming  that  the 
quasi-TEM  waves  are  non-dispersive,  and  each  mode  is  launched  and 
transmitted  without  a  distortion.  Assuming  that  the  starting  ampli¬ 
tude  at  (z*0)  is  denoted  by  a^Ct),  then  at  z  >  (4  the  amplitude 
becomes 

ai4(t)  *  ai0 

where  v^  is  the  velocity  of  the  ith  mode.  If  the  load  network  consists 
of  passive  elements  such  as  capacitor  and  inductors,  the  shape  of  the 
reflected  waves  will  be  distorted  as  compared  to  the  incident  waves. 


APPENDIX  (F) 


COUPLING  COEFFICIENTS  FOR  TRANSVERSE  ELECTROMAGNETIC  (TEM) 
MODE  EXCITED  BY  LOCALIZED  CURRENT  SOURCES 


Figure  (FI)  represents  an  electric  current  source  J  located  inside 
a  parallel  plate  waveguide  which  produces  outgoing  waves  that  carry 
energy  to  ports  Z^  and  Zj.  The  total  electric  and  magnetic  fields 
traveling  in  the  positive  z  direction  have  been  derived  by  Collin  [FI] 
as 


®  -i  S  z 

E  +  *  l  c+(e  +  e  )e  n 
n  n  nz 

n*l 


r  +  -  'jSn2 

H  ’  Ji  c“  0  +  “  8 


(FI) 

(F2) 


where  n  stands  for  the  nth  outward  propagating  mode,  c+  is  the  wave 

amplitude,  and  3 q  is  the  nth  mode  propagation  constant.  The  vectors 

e  and  h  are  normalized  transverse  electric  and  magnetic  fields  of  the 
n  n 

nth  quasi-TEM  mode  and  are  real  functions  normalized  to  unit  power  as 


s 


(F3) 


where  o  is  the  Kroneker  delta  and  s  is  the  cross  section  of  the  wave¬ 
guide.  Collin  has  used  the  reciprocity  theorem  for  Figure  (FI)  and 
has  obtained 


+ 

c 
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J  e 


j  3  z 
n 


dV 


(F4) 


where  the  volume  of  integration  encloses  the  current  J. 


Figure  (FI).  A  current  source  excitation  of  a  waveguide  [FI]. 


118 


The  above  method  can  be  used  In  aperture  problems  where  a  wire  is 
located  behind  a  conducting  screen  as  shown  in  Figures  (F2) .  The 
incident  field  from  exterior  region  produces  both  a  radiation  field 
and  a  guided  transverse  electromagnetic  wave  (TEM)  in  the  interior 

w 

region.  The  radiation  field  (Era{j)  radiates  into  a  free  space  and 
the  TEM  wave  propagates  along  the  wire.  In  order  to  use  the  method 
of  Collin,  one  must  assume  that  the  radiated  field  is  negligible. 
Kajfez  [F2]  has  attempted  to  justify  this  assumption  in  his  coupling 
formulation  [F2].  Davis  [F3]  has  also  compared  the  radiated  and  the 
transmitted  TEM  energy  for  such  problems  to  justify  the  assumption 
that  only  TEM  modes  need  to  be  considered. 

Shown  in  Figure  (F3)  are  the  electric  and  magnetic  dipoles  above 
a  closed  aperture  which  replace  the  aperture.  The  discussion  of  the 
dipole  equivalence  in  aperture  problems  is  formulated  in  Appendix  (B) . 
The  situation  in  Figure  (F3)  is  entirely  the  same  as  Figure  (FI)  with 
the  two  dipoles  and  Jge  being  the  effective  sources.  Now,  one 
can  apply  the  reciprocity  method  of  Collin  to  obtain 

c+-i  (h  •  Jm  -  e  •  J 8 )  e^  n  ds  (F5) 

n2'nsns; 

. 

Aperture 

where  the  integration  is  over  the  surface  of  the  aperture. 

The  computation  is  simplified  if  the  source  currents  are  assumed 
to  take  a  Dirac  delta  distribution  form  of  5(x)  and  5(z).  For  such  a 
distribution  equations  (Bll)  and  (B12)  become 
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ie 

s 

Substitution  of  equations 


"  cB  S(x  ~  xQ)  <$(z  - 
-  ce  6(x  -  Xq)  6(z  - 
(F6)  and  (F7)  in  (F5) 


Zo> 

(F6) 

zo) 

(F7) 

gives 

+  i  +3Vo 


cn  'I* 


h  (xft,  0)  c  -  e  (xrt,  0)  c 
_xn  0  mx  yn  0  eyj 


(F8) 


Similarly,  one  can  show  that 


1  /jenzo  r 


2  ® 


I  (x_ »  0)  c  -  a  (Xp.,  0)  c  i 
|_  xn  0  mx  yn  0  ey  J 


(F9) 


where  the  minus  on  h  _ (x_,  0)  and  in  the  exponent  are  due  to  the 

xn  (J 

reference  direction  for  propagation  in  the  (-z)  direction.  Equations 
(F8)  and  (F9)  completely  define  the  source  coupling  to  a  wire  behind 


an  aperture. 
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